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AN  MHD  INSTABILITY  PRIMER 


Introduction 

This  NRL  memo  report  provides  what  we  hope  is  a relatively  simple, 
self  contained  description  of  MHD  instabilities.  It  is  the  result  of  a I 

one  year  sabbatical  (Sept.  77  - Sept.  78)  one  of  us  (WMM)  spent  at  Culham 
Laboratory  in  Abingdon,  England.  The  principle  purpose  of  this  sabbatical 
was  to  learn  the  modern  theory  of  MHD  instabilities.  During  the  year 

J 

there,  the  two  of  us  worked  very  closely  in  our  attempt  to  learn  this 
field.  In  this  endeavor,  we  each  received  a tremendous  amount  of  help 
from  Dr.  John  Wesson,  also  of  Culham. 

Although  there  are  many  plasma  physicists  who  have  had  much  more 
experience  in  this  area  than  we  have,  we  feel  that  having  just  learned  it 
ourselves  makes  us  particularly  sympathetic  to  the  problems  someone  has 
in  trying  to  learn  about  MHD  instabilities  for  the  first  time.  Indeed  one 
of  the  problems  of  this  field  seems  to  be  that  it  is  extremely  difficult  t 
for  someone  to  learn  it  by  reading  the  primary  references.  This  has  led 
to  a situation  where  a relatively  small  number  of  researchers  have  gotten 
far  ahead  of  the  rest  of  the  plasma  physics  community.  Our  feeling  is 
that  scientific  endeavors  will  make  more  rapid  progress  if  knowledge  is 
broadly  based.  Thus  we  feel  that  there  is  a need  to  make  contact  between 
researchers  in  MHD  and  the  rest  of  the  community.  Our  intention  here  is 
to  help  provide  this  link  up;  in  doing  so  we  hope  to  provide  a service  not 
only  for  those  who  would  like  to  learn  this  area,  but  also  for  those  who 
already  know  it. 

One  deceptive  aspect  of  MHD  instabilities  is  that  the  simplest  ones 
are  extremely  easy  to  understand.  For  instance  the  instabilities  of  a Z 
pinch  to  sausage  and  kind  displacements  have  been  described  in  such  standard 
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texts  as  Jackson's  Classical  Electrodynamics.  However  more  complicated 
instabilities,  for  instance  in  a plasma  where  both  an  axial  and  azimuthal 
field  are  present  are  much  more  difficult  to  visualize;  but  they  are  also 
much  more  interesting  and  fun  to  learn. 

Although,  as  we  will  see,  there  is  a tremendous  variety  of  MHD 
instabilities,  all  those  which  we  will  study  (and  we  cover  most  of  the 
major  ones)  are  driven  by  one  of  two  mechanisms.  First  of  all  there  may 
be  a gravity  (or  something  equivalent  to  a gravitational  force)  which  is 
opposite  to  the  density  gradient.  This  can  drive  Raleigh-Taylor  type 
instabilities.  This  is  the  fundamental  driving  force  behind  both  ideal 
and  resistive  pressure  driven  modes  in  a reversed  field  pinch,  ballooning 
modes  in  a takamak,  rotationally  driven  modes  in  a 9 pinch  and  the  mirror 
instability  in  a magnetic  mirror  whose  field  lines  bulge  outward  away  from 
the  plasma.  Second,  a plasma  which  carries  a current  is  potentially 
unstable  because  current  elements  traveling  in  the  same  direction  attract 
each  other  and  would  like  to  all  clump  up  together.  Of  course  this  is 
not  so  simple  because  the  current  flows  through  a conducting  fluid  and, 
as  will  be  amply  demonstrated,  this  imposes  all  sorts  of  constraints. 
Nevertheless,  this  mutual  attraction  of  like  current  elements  is  the  basic 
mechanism  which  drives  free  surface  modes  in  a cylindrical  or  toroidal 
plasma;  tearing  modes  in  plane,  cylindrical  and  tokamak  geometry,  and 
the  internal  m = 1 kink  tearing  mode  in  cylindrical  or  tokamak  geometry. 

There  are  other  views  we  have  on  this  area  which  affected  the  material 
we  chose  to  cover.  First  of  all,  we  de-emphasized  both  the  energy  principle 
and  also  modes  in  a plasma  with  a free  surface.  At  least  for  one  dimensional 
configurations  (that  is  plasmas  where  all  variation  is  in  one  direction), 
which  we  emphasize,  it  is  really  no  simpler  to  utilize  the  energy  principle 
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than  it  is  to  solve  for  the  eigenfunction  and  eigenvalue.  This  is  especially 
true  now,  where  second  order  ordinary  differential  equations  can  be  so  easily 
solved  numerically.  However  we  could  hardly  write  this  memo  and  make  no 
mention  of  the  energy  principle  since  it  is  widely  used.  Therefore  one 
chapter  is  devoted  to  it. 

As  far  as  free  surfaces  are  concerned,  first  of  all  it  is  rare  that 
plasmas  have  free  surfaces;  usually  experimental  plots  of  say  density  or  J 

temperature  profile  usually  show  them  going  smoothely  to  zero  with  radius. 

Secondly,  free  surfaces  are  not  particularly  difficult  to  understand,  but 
applying  boundary  conditions  across  them  in  their  unperturbed  and  perturbed 
state  can  involve  a great  deal  of  mathematical  complexity.  Therefore,  they 
do  not  seem  to  be  worth  expending  a great  deal  of  effort  on  in  a manuscript 
like  ours  which  attempts  to  emphasize  physical  principle,  not  mathematical 
detail.  Hence  we  deal  with  free  surfaces  in  only  one  chapter  in  which  we 
drive  some  fundamental  stability  requirements  for  tokamak  plasmas. 

One  thing  which  we  attempt  to  emphasize  however  is  magnetic  reconnection. 

As  we  will  see  shortly,  in  ideal  MHD,  each  magnetic  field  line  maintains  its 
integrity  and  behaves  rather  like  a string  which  threads  the  fluid  and  cannot 
break.  However  there  are  flow  patterns  which  force  two  magnetic  field  lines 
together.  Often  non  ideal  effects,  for  instance  resistivity,  can  cause  the 
field  lines  to  break  and  reconnect.  Not  only  is  reconnection  very  much  in 
vogue  these  days,  it  is  also  extremely  interesting  in  its  own  right.  Even 
if  a plasma  physicist  does  not  want  to  work  on  tearing  modes,  it  might  still 
be  worth  learning  about  them  for  the  sheer  joy  of  it. 

The  outline  of  the  remainder  of  this  memo  is  as  follows.  Chapter  II 
introduces  MHD  and  Chapter  III  proves  the  energy  principle.  Chapter  IV 
discusses  free  surface  modes  in  a cylindrical.  Chapter  V looks  at  gravity 
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driven  modes  and  shows  how  magnetic  shear  can  stabilize  them.  Chapter  VI 
shows  how  the  presence  of  resistivity  can  destabilize  gravity  modes  which 
have  been  shear  stabilized.  Chapter  VII  examines  the  tearing  mode  in 
slab  geometry  and  shows  how  instabilities  which  are  forbidden  in  ideal 
MHD  are  possible  if  magnetic  reconnection  is  allowed.  Chapter  VIII  looks 
at  how  all  of  these  instabilities  transform  when  going  from  slab  to 
cylindrical  geometry,  and  Chapter  IX  briefly  discusses  the  additional 
complications  of  toroidal  geometry. 
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II.  Introduction  to  MHD 


The  Magnetohydrodynamic  equations,  which  form  the  subject  of  this  work 
are  taken  as 

9/ 


(mass  conservation). 


— ^ + v*  ^ y * o 


p Y + /°(M  v = ' 


+ f 


(II  1) 


(II  2) 


(momentum  conservation  assuming  scalar  pressure) 


y b = o 

6.) 

7*e= 

(<■) 

(II  3) 


(Maxwells  equations).  In  Maxwell's  equations,  the  displacement  current  is 
neglected;  instead  the  electric  field  is  related  to  the  current  through  Ohms 
law 


£ + 


(II  4) 


where  E + — x B is  the  electric  field  in  the  reference  frame  moving  with 

— c — 

the  (nonrelativistic)  fluid  velocity  and  n is  the  resistivity.  The  only 
other  quantity  needed  is  the  pressure.  We  will  assume  an  adiabadic  law 


Equations  (II  1)  through  (II  5)  constitute  a complete  description  of  the 
system.  The  unknowns  are  density,  pressure  and  the  three  components  of 
magnetic  field,  electric  field,  current  density  and  fluid  velocity,  fourteen 
in  all.  Equations  (II  1 and  5)  are  two  scalar  equations  for  p and  p. 
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Eqs.  (II  2,  3c  and  4)  are  three  vector  equations  for  the  three  components  of 
v,  J and  E,  and  finally  Eq.  (II  3 a and  b)  are  equations  for  the  solenoidal 
and  irrotational  parts  of  B. 

Usually,  it  is  convenient  to  eliminate  E and  J directly  by  using  Eqs. 

(II  3c  and  4).  Doing  so,  the  momentum  equation  and  Maxwell's  equation 
become 

h-  Yf  * iiP  &*?'>* 

(II  6) 


(II  7) 


Much  of  our  discussion  will  concern  perfectly  conducting  Fluids  (n  = 0)  so 
that  Eq.  (II  7)  is 


li 

2* 


(II  8) 


Equations  (II  1,  5,  6 and  7 or  8)  form  a complete  description  of  the  magnetized 
conducting  fluid.  Derivations  of  these  equations  as  wall  as  discussions  of 
their  validity  and  possible  extensions  (for  instance  including  tensor  pressure, 
thermal  conduction,  finite  Larmor  radius,  etc.)  have  been  discussed  in  many 
textbooks  on  plasma  physics.  We  simply  assume  these  equations  describe  the 
plasma  and  investigate  their  consequences;  particularly  we  focus  on  how  the 
magnetic  field  couples  to  the  fluid  motion. 

Before  doing  this,  it  is  worthwhile  to  quickly  review  what  magnetic 
field  lines  and  flux  tubes  are.  Any  vector  field  has  streamlines.  For 
the  magnetic  field,  these  are  the  solution  of 


jx  _ dz. 
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The  field  lines  are  then  everywhere  parallel  to  the  magnetic  field.  Let  us 
now  imagine  an  element  of  area  6A  which  is  parallel  to  a field  line  at  some 
point  s on  the  field  line.  If  the  magnetic  field  at  this  point  has  strength 
B(s),  the  flux  through  this  element  of  area  is  6ij>  = B(s)  6A  = B(s)«<5A.  One 
can  then  imagine  a tube  of  constant  flux  around  the  field  line.  If  the  flux 
is  d<f>,  the  area  of  the  flux  tube  as  a function  of  distance  along  the  field 
line  is  given  by 

A(s)~ 


The  amazing  thing  about  a flux  tube  is  that  in  a perfectly  conducting 
plasma  (n  = 0),  which  we  will  refer  to  as  ideal  MHD,  a flux  tube  is  convected 
with  the  flow.  It  is  a simple  matter  to  prove  this  from  Eq.  (II  7).  Inte- 
grating Eq.  (II  7)  over  a fixed  area  5A,  we  find 


J = / (els*  vjH).  3 


(ii  id 


where  f ds  is  an  integral  around  the  closed  line  bounding  5<j>.  However 
ds  x v dt  is  the  area  swept  out  by  an  element  ds  of  the  periphery  of  the 
flux  loop.  The  term  on  the  right  hand  side  of  Eq.  (II  11)  is  then  negative 
the  flux  ^?ept  out  by  the  circumference  in  its  trajectory.  Therefore,  the 
total  flux  through  a surface  area  moving  with  the  fluid  does  not  change 
time,  or 


JJ 


- o 


(II  12) 


Perhaps  a more  direct  way  to  see  this  is  to  note  that  if  q = 0,  Ohms  law 
Eq.  (II  4),  simply  says  that  E = 0 as  one  moves  with  the  fluid.  However 
in  the  inertial  frame  locally  moving  with  the  fluid.  Maxwells  equation  says 
that  the  rate  of  change  of  the  enclosed  flux  is  minus  the  loop  voltage, 
that  is  zero. 
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While  the  magnetic  lines  of  force  are  well  defined  by  Eq.  (II  9),  in 
general  there  is  no  unique  way  to  define  the  motion  of  lines  of  B in  a 
changing  medium.  However  in  a perfectly  conducting  fluid,  the  flux  is 

i 

frozen  into  the  flow.  Since  the  flux  tubes  can  be  regarded  as  bundles  of 
field  lines,  one  can  equally  look  upon  the  field  lines  themselves  as  being 
carried  along  with  the  flow.  We  will  now  examine  just  what  this  means. 

Since  V*B  = 0,  the  field  lines  have  no  start  or  finish,  but  either  close  j 

on  themselves  or  else  have  infinite  length.  Since  they  are  frozen  into  the 

flow,  the  field  lines  cannot  reconnect,  or  in  other  words,  the  topological 

properties  of  the  field  is  maintained.  That  is,  the  field  line  can  stretch 

and  bend,  but  it  cannot  change  its  topology. 

Let  us  illustrate  this  for  two  dimensional  motion.  Say  the  field  line 
initially  is  the  dotted  circle  shown  in  Fig.  (II  la).  Since  the  velocity 
field  is  a single  valued  function  of  r for  all  time,  there  is  no  way  that 
two  fluid  elements  initially  far  from  each  other  can  ever  occupy  the  same 
point;  to  do  so  would  mean  the  fluid  elements  pass  through  each  other, 
implying  a double  valued  velocity  field.  Therefore,  while  a complicated 
flow  pattern  can  greatly  contort  the  field  line,  to  for  instance  the  solid 
line  in  Fig.  (II  la),  the  field  line  can  never  cross  itself,  it  always  has 
a single  inside  and  outside.  Thus  the  topological  properties  of  the  field 
line  are  maintained.  One  possible  flow  pattern  could  distort  two  nearby 
field  lines  to  a pattern  shown  in  Fig.  (II  lb).  While  the  two  points  near 
A can  be  arbitrarially  close,  the  field  line  still  must  maintain  its  integrity. 

However  as  we  will  see  in  later  chapters,  the  presence  of  non  ideal  effects, 
for  instance  resistivity,  can  relax  the  topological  constraint  so  that  the 
field  line  can  break  and  reconnect,  forming  that  pattern  in  Fig.  II  lc. 

Clearly  the  topology  has  changed  from  a simple  closed  curve  topology  to  a 
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figure  eight.  That  is  there  are  two  inner  regions  instead  of  one  and  two 
sets  of  closed  field  lines. 

As  another  example  the  topological  constraints  on  the  field  line 
motion,  we  will  consider  an  example  in  three  dimensions.  Consider  two 
field  lines  which  are  initially  parallel  to  each  other  shown  in  Fig.  (II  2a). 

(In  order  to  distinguish  which  field  line  is  in  front  of  the  other,  each 
line  is  shown  with  finite  width  and  different  shading).  Imagine  a flow 

1 

pattern  which  distorts  the  field  line,  but  which  renders  this  field  line 
periodic  in  space  with  periodicity  length  L,  and  which  has  v = 0 along  the 
dark  field  line. 

Since  the  lighter  field  line  is  frozen  into  the  flow,  it  cannot  wind 
around  the  dark  field  line  between  Z = 0 and  Z =.  L.  For  instance,  if  the 
fluid  winds  around  the  inner  field  line  at  say  Z = but  does  not  move 
at  Z = 0 and  Z = L,  the  shaded  line  winds  around  the  dark  line  as  shown  in 
Fig.  (II  2b).  Between  Z = 0 and  Z = /£,  the  shaded  line  winds  around  the 
dark  line,  but  between  Z = and  Z = L it  unwinds.  The  total  winding 
number  of  the  shaded  line  around  the  dark,  between  0<Z<L,  is  preserved  at 
zero. 

However  another  type  of  fluid  motion  might  give  rise  to  the  field  lines 
shown  in  Fig.  (II  2c).  The  topology  of  the  shaded  line  is  the  same  as 
shown  in  Fig.  (II  2a  and  b) . That  is,  if  one  pulled  the  field  line  at  Z = 0 
and  Z = L,  it  would  snap  back  to  that  shape  in  Fig.  (II  2a).  Notice  though 
that  near  the  point  marked  A,  two  parts  of  the  field  line  are  forced  close 
together.  Again,  the  presence  of  non-ideal  effects  could  cause  the  field 
line  there  to  break  and  reconnect  as  shown  in  Fig.  II  2d  where  now  two  sets 
of  field  lines  are  produced,  the  main  field  line  which  now  loops  around  the 
axis,  and  an  additional  circle  which  also  loops  the  axis.  Thus,  in  ideal 
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MHD  if  a field  line  does  not  wind  around  another  initially,  it  never  does. 
However  if  non-ideal  effects  (for  instance  resistivity)  are  allowed,  the 
field  lines  may  break  and  reconnect  and  then  wind  around  another  field  line 
which  it  did  not  initially  encircle.  In  other  words,  new  magnetic  axes 
can  be  generated. 

To  summarize,  the  field  lines  in  a magnetized  fluid  can  be  regarded 
as  strings  which  thread  the  fluid  and  go  wherever  the  fluid  goes.  They  can 
stretch  and  bend  but  cannot  break  or  reconnect  in  ideal  MHD.  However  if 
non-ideal  effects  are  allowed,  the  field  lines  can  break  and  reconnect. 
Obviously,  however,  since  div*B  = 0,  a field  line  cannot  break  unless  it 
reconnects  instantaneously  with  another  part  of  the  field  line  (or  with  a 
different  field  line).  As  we  will  see,  there  are  types  of  fluid  motion 
which  tend  to  force  different  portions  of  field  lines  together,  as  shown 
at  points  A in  Fig.  (II  lb  and  II  2c).  In  this  case,  often  a very  small 
amount  of  resistivity  can  cause  reconnection  at  these  points.  In  other 
words,  given  a choice  between  evolving  toward  a very  complicated  structure 
with  the  same  topology,  or  a simple  structure  with  different  topology,  a 
field  line  in  a real  plasma  will  often  choose  the  latter. 

We  now  turn  to  a study  of  the  magnetic  forces  exerted  on  the  plasma. 
The  magnetic  force  per  unit  volume  is  given  by 


Fr- 


(II  13) 


which  is  clearly  always  perpendicular  to  both  B and  J.  A more  convenient 
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In  Eq.  (II  14b)  above  V means  the  portion  of  the  gradient  which  is 
perpendicular  to  B and  i,  is  a unit  vector  in  the  direction  of  B. 

The  two  terms  in  Eq.  (II  14b)  have  simple  interpretation.  The  first 
term  shows  that  the  magnitude  of  B2  acts  like  a pressure  in  a direction 
perpendicular  to  B.  That  is  a gradient  in  B2  exerts  a force  which  pushes 
the  plasma  toward  regions  of  lower  B2.  To  interpret  the  second  term,  note 
that  (i,  *V)i,  = -R/  where  R is  the  radius  of  curvature  of  the  field  line 
and  i^  is  a unit  vector  pointing  toward  the  center  of  curvature.  Thus  if 
the  field  lines  are  bent,  there  is  a force  exerted  on  the  plasma  of 

g2 

magnitude  and  directed  toward  the  center  of  curvature.  This  latter 

force  acts  rather  as  if  the  field  line  were  a rubber  band.  If  a stretched 
rubber  band  is  bent,  a force  is  exerted  which  tends  to  snap  it  back  to  a 

straight  line.  This  then  corresponds  to  the  way  a magnetic  field  line  acts. 

The  main  difference  is,  of  course,  that  a stretched  bent  rubber  band  exerts 
a strong  force  along  its  length  as  well  as  perpendicular  to  itself;  the  bent 

magnetic  field  line  only  exerts  this  force  perpendicular  to  itself  but  not 

along  its  length.  Fig.  (II  3)  illustrates  the  force  exerted  on  the  plasma 
by  the  two  terms  in  Eq.  (II  14b). 

Clearly,  before  we  can  concern  ourselves  with  MHD  instabilities,  a first 

step  is  an  examination  of  MUD  equilibria.  In  equilibrium,  with  V = — = 0, 

0 t 

the  momentum  conservation  equation  is 


-L  Ty'B  = - 
c - — 


-J-V' + 
s<  - 


yir  v - ' — 


(II  15) 


There  are  several  immediate  consequences  of  Eq.  (II  15).  First  of  all,  it 
is  easy  to  show  that  the  pressure  must  be  constant  along  a field  line; 
simply  take  dot  product  of  Eq.  (II  15)  with  i^  and  get  the  result 
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(II  16) 


(Vs)  P =0 

so  the  pressure  is  constant  along  a field  line.  Therefore  any  open  ended 

device,  for  instance  a magnetic  mirror,  where  field  lines  end  on  walls,  cannot  be  in 

MHD  equilibrium,  at  least  for  scalar  pressure.  It  is  for  this  reason  that 

the  research  on  open  ended  devises  usually  concerns  itself  principly  with 

Vlasov  equilibria  and  stability  or  else  with  tensor  pressure. 

For  the  remainder  of  this  chapter  we  will  concern  ourselves  with 
devices  with  enclosed  field  lines,  for  instance  tokamaks,  reversed  field 
pinches,  or  infinitely  long  cylinders.  As  a field  line  goes  around  a 
toroidal  machine,  there  are  three  possibilities;  first,  it  may  close  on 
itself  after  one  or  more  transits;  second,  it  may  trace  out  a surface,  and 
third,  it  may  ergodigally  fill  a volume,  which  may  be  either  the  entire 
volume  of  the  torus  or  a portion  of  it.  In  MHD  equilibrium,  the  pressure 
is  constant  along  each  magnetic  line,  surface  or  volume,  whatever  the  case 
may  be.  Clearly,  if  the  field  lines  fill  the  volume  of  the  device,  no  MHD 
equilibrium  is  possible. 

In  the  case  of  reverse  field  pinches  or  tokamaks,  the  field  lines  form 
surfaces.  From  Eq.  (II  15),  it  is  possible  to  show  that  the  current  lines 
are  also  in  the  flux  surfaces.  To  do  so,  take  the  dot  product  with  J which 
shows  that  P = constant  along  the  lines  of  J.  However  if  P varies  from  one 
flux  surface  to  the  next,  then  lines  of  J must  also  lie  in  the  surfaces  of 
constant  ? that  is  the  flux  surfaces. 

In  a general  toroidal  equilibrium,  it  is  usually  convenient  to  choose 
as  one  co-ordinate  surface, the  flux  surface.  This  gives  rise  to  the  Grad- 

Shafranov  equation  and  is  discussed  in  Appendix  . For  our  purposes 

here,  we  only  write  the  pressure  balance  equation  in  cylindrical  geometry 
with  variation  only  in  r.  These  are 
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r 


0-  7~ 


I 


4* 


(II  17) 


Once  p and  are  given,  Eq.  (II  17)  can  be  solved  for  the  appropriate  BQ 
for  equilibrium. 

Having  discussed  the  steady  state  equilibrium,  we  now  proceed  to  a 
discussion  of  the  dynamics  by  examining  what  sort  of  wave  motion  is  allowed 
in  a uniform*  magnetized,  current  free  conducting  fluid.  Denoting  an 
unperturbed  (equilibrium)  quantity  with  a subscript  zero  and  a perturbed 
quantity  with  no  subscript,  the  linearized  equations  of  motion  are 

'f*r 

-f  o„  v-y  -o  0 (ii  is) 

* v x v x 5e 


where  we  have  made  use  of  the  adiabatic  relation  between  p and  p.  Taking 


— of  Eq.  (II  18)  and  substituting  from  Eq.  (II  18  b and  c),  we  find  a 

O t 

single  vector  equation  for  V.  Assuming  V *»  ik  and  = - im,  it  is 

= -<■/>(>•*) + v;  (,*[£’ [**(***)}] 


(II  19) 


where  c.g  is  the  sound  speed,  cg 


B 


0/\ 


('V. )‘ 


/ 2 


and  is  the  Alfven  speed 


Also  B_  is  assumed  to  be  in  the  z direction.  Since  the 


'A  / \ 4irpQ  ' -o 

plasma  is  isotropic  in  the  x-y  plane,  we  may  take  k to  be  in  the  yz  plane 
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without  loss  of  generality.  Then  taking  the  component  forms  of  Eq.  (II  19) 


gives 

U)VX  - 


U)  vz  = 


(0 

(II  20) 


As  is  clear  from  Eq.  (II  20a),  the  x velocity  decouples  and  this  is 
the  shear  Alfven  wave  with  dispersion  relation 


c&l  9 


(II  21) 


k being  the  magnitude  of  k and  9 being  the  angle  between  B and  k.  The  other 
modes  involve  the  coupled  y and  z motion.  A straightforward  calculation  from 
Eqs.  (II  20  b and  c)  gives  the  result 


(II  22) 


For  9 » 0,  the  two  roots  are  w2/k2  “ c2  and  V2,  while  for  9 « /£,  the  two 

c 2 

roots  are  Cg  + v \ a"d  oi2/k2  * c|  cos29/(l  + a/j2^  **  0.  A polar  plot 

A 

/kV^  for  the  three  roots  is  shown  in  Fig.  (II  4)  for  the  case  of  VA  » 3cg. 
Figure  (II  4)  shows  that  waves  in  a magnetized  fluid  have  some  interesting 
properties.  Since  the  shear  Alfven  wave  (Eq.  (II  21))  has  V * V * 0,  it 

y z 

does  not  compress  the  plasma.  At  0 a 0,  the  sound  wave  does  compress  the 
plasma,  but  by  following  this  root  around  to  9 - % where  u = 0,  one  can  show 


that  there  is  no  plasma  compression  there.  On  the  other  hand,  at  9 = 0, 


the  wave  which  has  ^ does  not  compress  the  plasma;  however  following 

the  root  around  to  9 = %,  where  = Va  + Cs’  we  See  that  the  plasma  is 

compressed.  Thus  if  V > c , there  is  no  such  thing  as  a pure  sound  wave. 

What  starts  out  as  a compressional  wave  at  9 = 0 ends  up  as  a shear  wave  at 

9 = /C,  and  visa  versa.  However  if  c > V.,  the  sound  wave  is  compressional 
Z S A 

at  all  angles. 

We  now  give  simple  physical  pictures  first  for  the  shear  Alfven  wave 
for  9=0  and  then  for  the  compressional  Alfven  (magnetosonic)  wave  at 
9 = /^.  Imagine  that  each  fluid  element  is  displaced  in  the  x direction 
an  amount 

(II  23) 


Since  the  field  line  is  frozen  into  the  flow,  Eq.  (II  23)  above  is  also  the 
equation  for  a field  line.  The  perturbed  field  B i is  perpendicular  to 

X“*X 

the  equilibrium  field,  BQiz,  so  there  is  no  first  order  change  in  the 
magnitude  of  B.  Hence  the  only  force  on  the  fluid  arises  from  the  bending 
of  the  field  line,  as  discussed  after  Eq.  (II  14).  If  terms  of  order  xq2 
are  neglected,  the  reciprocal  of  the  radius  of  curvature  of  the  field  line 
is 


ix 


(II  24) 


Therefore  the  force  per  unit  volume  on  each  fluid  element  is 


a7-  . 

--  iX  Zo 


(II  25) 


which  is  just  minus  a constant  times  the  displacement  of  the  fluid  element. 
That  is,  each  fluid  element  performs  simple  harmonic  oscillation  with 
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frequency  given  by  to2  = k2BQ2/4irpo  perpendicular  to  the  equilibrium 
magnetic  field.  The  phase  speed  of  this  wave  is  the  Alfven  speed. 

One  can  picture  this  oscillation  in  terms  of  the  magnetic  field  lines 
constituting  a series  of  strings  which  permeate  the  plasma.  Imagine 
pulling  on  these  strings  with  a tension  of  B2^  per  unit  area.  Then  if 
the  strings  are  bent,  they  will  tend  to  snap  back  and  oscillate  about 
their  equilibrium  position  with  frequency  k2B2/47rp, 

Therefore  as  far  forces 

perpendicular  to  themselves  are  concerned,  the  magnetic  field  lines  are 
like  strings  with  tension  B2^  per  unit  area.  They  are  different  however, 
in  that  unlike  a stretched  string  there  is  no  force  along  a magnetic  field 
line. 


We  now  turn  to  an  examination  of  wave  motion  across  a magnetic  field. 
Say  that  Bq  is  in  the  z direction  and  the  displacement  and  wave  number  are 
both  in  the  y direction.  In  this  case,  the  field  lines  remain  straight 
and  do  not  bend.  Therefore,  according  to  Eq.  (II  14b),  the  force  density 
on  the  plasma  is  minus  the  gradient  of  the  scalar  magnetic  pressure.  This 
perturbed  magnetic  pressure  is 


U 


'do 


(II  26) 


where  we  have  used  the  fact  that  the  field  is  frozen  into  the  plasma  so 
p/  = (Also  this  can  be  easily  verified  from  Eqs.  (II  18b  and  c).). 

P * D 

0 zo 

Hence  the  magnetic  scalar  pressure  acts  like  a fluid  with  y = 2.  The 
perturbed  fluid  pressure  is  of  course  ^p  so  that  the  perturbed  total 

pressure  in  the  plasma  is  given  by 


' To-r 


(12  + 

v r- 


(II  27) 
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Now  it  is  clear  that  perpendicular  to  the  field,  waves  propagate  like  sound 
waves,  except  that  the  magnetic  field  adds  an  extra  'springiness'  to  the 


plasma. 

We  conclude  this  chapter  with  a discussion  of  what  sort  of  plasma 
motion  an  MHD  instability  is  likely  to  generate.  The  modes  we  have  discussed 
all  have  u)2  > 0.  However  for  instability,  obviously  w2  < 0.  Now  imagine 
that  there  is  some  physical  effect  which  perturbs  the  system  so  as  to  drive 
the  plasma  toward  instability.  Clearly  this  effect  will  manifest  itself  by 
lowering  a)2.  However  if  a>2  is  large  to  begin  with,  instability,  in  general, 
will  not  result,  but  instead  just  a lowering  of  the  frequency.  Clearly,  the 
most  likely  place  for  instability  is  where  w2  = 0,  or  for  incompressible 
flow  with  k*B  = 0,  according  to  Fig.  II  4 and  the  discussion  following 
Eq.  (II  22).  In  other  words,  if  k*B  # 0,  the  unstable  flow  will  couple  to 
shear  Alfven  waves,  which  is  a stabilizing  effect;  and  if  7»V  f4  0,  it  will 
couple  to  sound  or  magnetosonic  waves,  which  is  also  stabilizing.  In  most 
of  the  rest  of  this  book,  we  focus  on  perturbed  plasma  motion  which  is 
incompressible.  This  does  not  mean,  of  course,  that  the  plasma  is  impossible 
to  compress,  but  rather  that  the  motions  we  consider  have  phase  speed  much 
less  than  the  magnetosonic  speed  and  are  therefore  decoupled  from  compressional 
motion.  This  is  no  different  from  any  other  fluid  which  is  assumed  to  be 
incompressible.  For  instance  with  regard  to  paddling  a canoe,  water  is 
incompressible;  however  sound  does  propagate  through,  and  does  compress  it. 
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Fig.  II  4 — A polar  plot  of  phase  velocity  versus  angle  to  the  magnetic  field 
for  the  three  branches  of  MHD  oscillation. 
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III.  The  Energy  Principle  | 

A complete  solution  of  the  linear  stability  of  a plasma  for  a given  1 

magnetic  field  configuration  requires  the  determination  of  eigen-functions 
and  eigen-values  of  the  linearized  equations.  For  all  but  the  simplest 
geometries  this  proves  to  be  a problem  of  considerable  complexity.  It  is 

therefore  desirable  to  have  a procedure  for  deciding  the  question  of  stability  1 

/ 

which  does  not  require  the  determination  of  the  eigen-values  (or  characteristic 
frequencies).  For  the  problem  of  the  MHD  stability  of  a perfectly  conducting 
plasma  in  the  absence  of  an  equilibrium  flow,  the  Energy  Principle  provides 
just  such  a method. 

The  original  derivation  of  the  Energy  Principle  by  Bernstein  et  al 
assumed  that  the  eigen-functions  of  the  linearized  equations  formed  a complete 
set.  However,  it  was  pointed  out  by  Laval  et  al  (G.  Laval,  C.  Mercier,  and 
R.  Pellat,  Nucl  Fusion  _5,  156  (1965))  that  this  assumption  is  not  always  valid 
so  that  it  became  necessary  to  find  a proof  of  the  Energy  Principle  which 
did  not  rely  on  the  completeness  property.  Laval  et  al,  in  a very  elegant 
analysis,  provided  this  proof  assuming  only  that  the  linear  operator  was 
self-adjoint.  Here  we  shall  give  an  alternative  proof  of  the  Energy  Principle 
which  assumes  neither  completeness  nor  self-adjointness  but  instead  we  shall 
demonstrate  the  conservation  of  small  signal  energy  directly  from  the 
linearized  ideal  MHD  equations. 

Before  giving  this  proof  let  us  say  a few  words  about  the  small  signal 
energy  mentioned  above.  The  small  signal  energy  is  defined  entirely  in 
terms  of  the  fields  of  the  linear  theory  and  is  not  the  same  as  the  physical 
energy.  It  was  Sturrock  who  first  drew  attention  to  the  significance  of  the 
small  signal  energy  (or  "pseudo-energy"  as  he  called  it)  for  the  problem  of 
linear  stability.  It  is  not  obvious  that  this  small  signal  energy  will  be 
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conserved  in  a conservative  system  simply  because  the  physical  energy  is 
conserved  since  the  exact  fields  have  been  expanded  in  a series  and  truncated 
at  the  first  order. 

We  shall  now  derive  a generalization  of  the  Poynting  theorem  from  the 
linearized  MHD  equations  and  from  this  obtain  the  conservation  of  small 

i 

signal  energy.  The  proof  of  the  Energy  Principle  is  then  completed  using 

the  arguments  given  by  Laval  et  al.  ^ 

Let  us  consider  a plasma  in  which  the  pressure  is  isotropic  and  which 

is  bounded  by  a rigid,  perfectly  conducting  wall,  where  the  boundary  conditions 

are  the  following 

» « * 3B 

n-V  = 0;  n X E =0;  n - =0 

where  g is  the  unit  normal  to  the  boundary.  The  analysis  is  easily  generalized 
to  other  boundary  conditions  (e.g.  plasma-vacuum  boundary).  The  linearized 
equations  of  the  ideal  MHD  model  are 


ay 


Vp  + f r ^ * ? 


(hi.  i) 


- £> 


(HI.  2) 


(in.  3) 
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(III.  4) 


^ -*■  (v.  ?)  P 


V*f  - ' 


J- 

c 


3 8 

3* 


(III.  5) 


V*S 


(III.  6) 


(in.  7) 


where  fields  with  subscript  zero  are  equilibrium  quantities  and  the  linearized 
variables  are  written  without  subscripts.  As  already  mentioned  there  is  no 
equilibrium  flow  of  the  plasma.  We  now  scalar  multiply  equation  (III.  1)  by 
v to  obtain 


v*  Yr 


(III.  8) 


using  the  relation 

V-  (pv)  - p <7-  v + v.  7 p 

equation  (III.  8)  can  be  written  in  the  form 


3<f 


- - v-(pyj  +•  1.^5,)+  /(Sr.*  ?J 


(III.  9) 


Next,  scalar  multiple  equation  (III.  6)  by  E, 


equation  (III.  5) 


by 


cB 

4tt 
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to  obtain 


7-  E-V  x B - 7-  B*V  x E 
4tt  - - - 4tt  - 


h <S-*j  + J-s 
3t  8^ 


This  can  be  written  as 


B-  ? 
8*  . 


+ V 


V ^ * 0 

V ’7/' 


(III.  10) 


Substituting  for  g from  equation  (III.  3)  into  the  term  J • E equation  (III.  10) 
becomes 

c e / B 

I - — V * I <-J  r L-Jo  ) 

(III.  11) 


uiuco 


where  we  have  made  use  of  the  vector  identy  J (vx  Bq)  « ■' *(B  x J).  Adding 
equations  (III.  9)  and  (III.  11),  we  obtain 


(III.  12) 


r* &***£)  + + r*) 

= p ?-y  * y-  ( J". 

Now  consider  the  term  p7*v.  With  the  aid  of  equations  (III.  2)  and  (III.  4), 
we  f ind 

£‘-v  ° - 

Substituting  equation  (III.  13)  into  (III,  12)  we  have 


(III.  13) 


X ) v.v  J- 

rf°~  ~ * ** 


(III.  14) 


s v • ( ^ * h)  - _X-  (y  • y ) p* 

In  order  to  obtain  the  final  form  of  the  conservation  equation  we  must 
introduce  the  linear  displacement  vector  £ defined  by 

XL  - v 

atf  ' - cm.  15) 
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This  enables  us  to  integrate  equation  (III.  5)  in  time  to  obtain  B in  terms  of 
§ thus 

B ' ffo)  (III.  16) 


where  we  have,  of  course,  made  use  of  equation  (III.  3).  We  can  also  obtain 
p in  terms  of  £ by  integrating  equation  (III.  13)  in  time  to  obtain 


P*  ' (?’*)?.- 


(III.  17) 


Now  consider  the  two  terms  on  the  right  hand  side  of  equation  (III.  14).  First, 

y- s.-cU'?) 

Since  2 depends  linearly  on  £ and  Jq  is  independent  of  time  we  may  write 


Sisf'?)--  3(* 


so  that 


(III.  18) 


Next,  consider  the  second  term 


Z.  If. 


• * f* 


Since  p depends  linearly  on  5 through  equation  (III.  17)  and  pQ  is  independent 


of  time  we  may  again  write 


(y-s) 


37* 


P°"  -Z*  l * Vpo-  - 1 ' \ (III.  19) 

We  may  now  substitute  (III.  18)  and  (III.  19)  into  equation  (III.  14)  to  obtain 
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i _ 

a. 


y 


£~  * B 


+■ 


(III.  20) 


This  is  the  required  energy  conservation  relation  and  is  the  generalized 
Poynting  theorem  for  linearized  ideal  MHD.  The  term  inside  the  first  set 
of  parentheses  is  the  total  energy  density  of  the  perturbation  per  unit 
volume  and  the  second  set  of  parentheses  represents  the  energy  flow  out  of 
this  unit  volume.  The  energy  flow  terms  will  be  easily  recognized  as  the 
Poynting  vector  and  the  convection  of  the  plasma  internal  energy  by  the 
perturbed  plasma  motion. 

Equation  (III.  20)  is  not  yet  in  its  most  familiar  form.  Let  us  combine 


the  pair  of  terms  containing  the  perturbed  pressure 

1 JL  3"  fr.  (?* 


* fo 


I ' 3 p.  - 


Substituting  for  p from  equation  (III.  17)  we  obtain 

* T.  >■  V £ 


-L  JL.  J-  -2-  f • V 

ypo 


a * yp,  - ’P°  = ) 


(III.  21) 


Substituting  (III.  21)  info  (III.  20)  we  obtain  the  final  form  for  the  small 
signal  energy  conservation  equation 

2-  \ j-Ay-  y * + £(f 

jt'  [ 1 

+ \I ' * F ~ J ' O (in.  22) 
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We  now  integrate  this  aquation  over  the  whole  plasma  to  obtain 

£{!*/**'**  v-ijd’r 


£ / B * 


Pyj-JS 


= o 


where  dS  = nS  and  S is  the  surface  bounding  the  volume  of  integration. 

Since  we  have  assumed  a perfectly  conducting,  rigid  wall  in  contact  with 
the  plasma  the  boundary  conditions  ensure  that  the  surface  integral  vanishes. 
We  are  then  left  with  the  result  that 


A ( K*  <fw)  = O 

^ (III.  23) 


where  K = J pQV*vdV  is  the  total  kinetic  energy  of  the  plasma  5W  is 
immediately  identified  as  the  potential  energy  of  the  perturbation  and  is 
given  by 


The  form  of  6W  given  in  equation  (III.  24)  is  in  agreement  with  that  first 
given  by  Bernstein  et  al. 

We  can  already  draw  a number  of  conclusions  at  this  stage  of  the 
analysis  with  the  aid  of  the  equation  for  the  conservation  of  energy.  Since 
we  assumed  zero  equilibrium  flow,  it  follows  that  K is  always  positive 
definite.  Taking  as  a definition  of  instability  an  unbounded  increase  of  K 
in  time,  it  is  clear  from  equation  (III.  23)  that  if  6W  > 0 the  system 
must  be  absolutely  stable.  It  is  also  clear  from  equation  (III.  23)  that 
in  order  to  have  instability  6W  < 0 such  that  | 6W | grows  in  time  so  as  to 
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balance  exactly  the  increase  in  K.  To  complete  the  proof  of  the  Energy 
Principle  we  now  follow  the  argument  given  by  Laval  et  al  in  order  to  show 


that  if  a 5 and  B = V x £ x Bq  can  be  found  such  that  6W  < 0 the  plasma 
is  unstable. 

Equation  (III.  1)  can  be  written  as 


(III.  25) 


where 


F(S)  = V ( ?•  v ?■>)+  * (*F‘ 

+ 

The  next  step  is  to  define  the  virial  1(5)  defined  by 

Ilf)--  jrfptSSr 

Differentiating  I twice  with  respect  to  time 

i-  ffai-h 


(III.  26) 


(III.  27) 


(III.  28) 


The  first  term  on  the  right  hand  side  is  twice  the  kinetic  energy  and  the 
second  term  can  be  written  in  terms  of  F(5)  with  the  aid  of  equation  (III.  25) 
giving 


j * °?k  i-  f i * 


(III.  29) 


With  the  aid  of  some  vector  algebra  it  is  straightforward  to  show  that 


J S-r(f)d'r*  - jsw 

where  6W  is  the  potential  energy  defined  in  equation  (III.  24). 
form  of  equation  (III.  29)  can  now  be  written 


(III.  30) 

The  final 


J = JK  - <2  SW 


(III.  31) 
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We  now  assume  that  there  is  some  displacement  n such  that  6W(n)  < 0.  We 
write  this  explicitly  as 


SV(y)'-  ^ 2 -(?) 


(III.  32) 


where  w > 0 and  I(n),  the  virial,  is,  of  course,  positive  definite.  The 
displacement  vector  C>  which  is  a solution  of  equation  (III.  25),  is  now 
chosen  to  satisfy  the  following  initial  conditions 

£ - 7 j i ~ ^7  (hi.  33) 

Since  the  total  energy  Q of  the  perturbation  must  be  a constant  we  may 
calculate  its  value  with  the  aid  of  equations  (III.  23).  Thus 


Q = K + 6W 

= 1(0  + <SW(§) 


= u)2I(ri)  + 6W(n) 

so  that  Q = 0.  Using  this  result  we  may  now  eliminate  6W  from  equation 
(III.  31)  to  obtain 


(III.  34) 


We  now  relate  I to  I and  I,  By  the  definition  of  I, 

• « 

I 


l f f>.  I 


Using  Schwarz's  inequality,  we  may  write 


(?  # (/>•  f-ft'O  ([  p-  J3-r) 


and  therefore  I2  < 4IK.  With  the  aid  of  equation  (III.  34)  the  inequality 
becomes 

* J -r  Y 

I K X (III.  35) 
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Now  I > 0 at  t = 0.  Integrating  equation  (III.  34)  we  obtain  I = I + f C4Kdt 

o Jq 

where  the  subscript  zero  indicates  that  the  quantity  has  been  evaluated  at 
t 

Thus 


0.  We  may  evaluate  the  constant  I with  the  aid  of  the  initial  conditions. 


r oPu>  T0  >0 


(III.  36) 


We  therefore  find  that  I > 0 so  that  I > 0 for  t > 0.  Since  II  > 0 for 

o 

t > 0 we  do  not  alter  the  sense  of  inequality  (III.  35)  by  dividing  through- 


out by  this  quantity,  so  that 

I . J 


(III.  37) 


Integrating  for  t > 0 we  have 

X,  % (■  JL  'x/Xo 

Using  inequality  (III.  36)  we  obtain 
-£•  > cPui 

I 

Again  integrating  for  t > 0 we  have 

1 ^ To  ^ 


(III.  38) 


(III.  39) 


(III.  40) 


Thus  I grows  at  least  as  fast  as  exp(2wt)  and  since  I is  quadratic  in  £ it 
follows  that  the  displacement  E,  will  grow  at  least  as  fast  as  exp(wt).  This 
then  completes  the  proof  that  for  any  displacement  which  makes  6W  < 0 
instability  will  always  occur. 

Let  us  conclude  this  chapter  by  considering,  in  very  simple  terms, 
the  physical  significance  of  the  potential  energy  <$W  defined  by  equation 
(III.  24).  We  have  seen  that  the  plasma  will  be  stable  when  6W  > 0.  The 
expression  for  6W  then  shows  that  the  energy  required  to  perturb  the  equilibrium 


A 
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magnetic  field  (either  field  line  bending  or  compression)  is  positive  definite 

and  therefore  stabilizing.  The  third  term  is  the  expression  for  6W  is  also 

positive  definite  so  that  incompressible  perturbations  may  be  expected  to  be 

the  most  unstable.  The  second  and  fourth  terms  in  the  expression  for  6W  are 

the  potentially  destabilizing  ones.  The  driving  mechanism  in  the  first  of 

these  being  due  to  the  equilibrium  current  and  is  the  second  to  the  equilibrium 

pressure  gradient.  To  apply  the  energy  principle  one  can  do  one  of  two  things. 

First  of  all,  one  can  try  various  different  displacements  and  see  whether  the 

6W  can  be  made  negative.  For  instance  the  well  known  stability  requirement 

f di. 

that  / pg  must  decrease  as  the  pressure  decreases  can  be  derived  by  examining 
the  effect  of  a displacement  which  interchanges  two  adjacent  flux  tubes.  Of 
course  this  only  gives  a necessary  condition  for  stability;  there  may  be 
other  untried  displacements  which  give  instability.  Alternatively  one  could 
attempt  to  minimize  6W  by  applying  the  variational  principle.  If  the  minimum 
value  of  6W  > 0,  for  a displacement  which  for  instance  has  J |s|2d3r  = l, 
then  the  plasma  is  stable,  and  visa  versa.  Thus  one  can  derive  a necessary 
and  sufficient  condition  for  stability.  This  is  the  approach  used  by  Newcomb 
(W.A.  Newcomb,  Annals  of  Physics,  10,  232  (I960))  in  analyzing  the  stability 
of  a diffuse  linear  pinch.  For  an  illustration  of  the  Energy  Principle  in  a 
number  of  simple  situations  the  interested  reader  may  refer  to  the  review  of 
hydromagnetic  stability  by  Kadomtsev.  (B.  B.  Kadomtsev,  Rev.  Plasma  Physics, 

Vol.  2,  Edited  by  M.  A.  Leontovich,  Consultants  Bureau,  New  York  (1966)  p.  153). 
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IV.  Free  Surface  Modes  in  a Cylindrical  Plasma 


A very  important  class  of  unstable  modes  occur  when  the  plasma  does  not 
extend  up  to  the  wall  of  the  metal  chamber.  In  this  case,  the  plasma  has  a 
free  surface  which  can  undergo  unstable  helical  (or  kink)  perturbations.  In 
this  chapter,  we  work  out  the  properties  of  two  types  of  free  surface  modes, 
first  where  the  current  is  carried  on  the  surface,  and  second  where  the 
volume  current  density  J ^ is  uniform  inside  the  plasma.  The  former  instability 
is  responsible  for  the  famous  Kruskal-Shafranov  requirement  for  . gross  stability 
q(a)>l  where  a is  the  radius  of  the  free  surface.  The  later  instability  is 
responsible  for  unstable  displacements  of  the  free  surface  wherever  q(a)<I5j1 
where  m and  n are  integers.  Also  q(r)  - rB  /RB  , R being  the  major  radius. 

We  begin  with  the  former  type  of  instability,  that  is  where  all  the  current 
flows  along  the  free  surface.  Let  us  now  define  our  model  more  precisely. 

We  shall  consider  an  infinitely  long  cylindrical  plasma  whose  axis  is  taken 
to  be  the  z direction.  The  plasma  extends  radially  from  r = 0 to  r = a, 
beyond  which  there  is  a vacuum.  For  simplicity,  assume  the  vacuum  region 
extends  to  infinity  and  that  the  plasma  surface  current  flows  only  in  the 
z direction.  Then,  Maxwell's  current  equation  gives  the  following  boundary 


conditions  for  equilibrium  magnetic  fields 


(IV  1) 


Vtr 


C 


(IV  2) 


where  Jqz  is  the  axial  surface  current  and  the  superscripts  V and  p denote 

vacuum  and  plasma  fields  respectively.  Assuming  that  B = — = — = 0 

r do  aZ 

for  the  equilibrium  we  obtain  the  following  equilibrium  fields  from  Maxwell's 
equations 


B 


p 

o© 


-O 


(IV. 3) 


(IV. A) 


V-rra.  -Toz 
r c 


(IV. 5) 


1 


Since  the  magnetic  field  is  uniform  in  the  plasma,  then  so  must  be  the 
pressure.  Then  equilibrium  pressure  balance  at  the  plasma  vacuum  interface 
is  given  by 


p.  = 


«?fT  -To1 


(IV. 6) 


Let  us  now  consider  perturbations  to  this  equilibrium,  which  are  assumed 

to  vary  as  f(r)  exp^i  (kz  + m9)  + yt^  (To  make  the  transition  from  cylindrical 

to  toroidal  systems,  we  simply  quantize  k by  taking  k ■ ^ where  R is  the 

major  radius  of  the  torus).  Here  and  henceforth  we  take  B , J and 

J 9o,  zo  zo 

m>0.  However  the  sign  of  k is  arbitrary.  In  general  there  is  only 
instability  for  k <0.  The  linearized  equation  of  motion  of  the  plasma  can  be 
written 


(IV. 7) 
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B *B 

where  pQ  is  assumed  to  be  uniform  and  P = P + ^ . In  obtaining  (IV.  7) 

we  have  used  the  fact  that  Bq  is  uniform  in  the  plasma.  Taking  the  divergence 


of  Eq.  (IV.  7)  we  obtain 


v)  * 


7*  B - V 


(IV.  8) 


We  now  assume  that  the  motion  is  incompressible  V*V  = V*B  = 0 so  Eq.  (IV.  8) 
reduces  to  the  simple  form 


v P 


(IV.  9) 


The  solution  of  this  equation  which  is  bounded  at  the  origin  is 


A T~C  AO 


(IV.  10) 


where  A is  an  arbitrary  constant  and  I is  the  modified  Bessel  function  of 

m 

the  first  kind. 

We  now  also  need  to  relate  the  perturbed  velocity  V to  P.  The  equation 
for  the  perturbed  magnetic  field  is 


jr  8-  V-x/vx  "B.  ) 


(IV.  11) 


Using  once  more  the  incompressibility condition  for  Bq  and  Maxwell's  equations, 
the  above  equation  reduces  to 


(IV.  12) 


Substituting  (IV.  12)  into  Eq.  (IV.  7)  we  obtain  V in  terms  of  P 


V * 


r Y ? 


(IV.  13) 
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where  C2  = B2  /4rrp  , the  Alfven  speed. 

A oz  o’ 

V V 

The  perturbed  field  in  the  vacuum  must  satisfy  V*B  = V x B = 0. 
v 

Putting  B = Vip,  then  it  satisfies  the  equation 
V2ip  = 0 


so 


< P = C Km(kr) 


(IV.  14) 


where  C is  a constant  and  Km  is  an  arbitrary  Bessel  function  of  the  second 
kind.  The  final  step  in  the  analysis  is  to  match  the  perturbation  in  the 
plasma  and  vacuum.  Since  we  have  two  arbitrary  constants  in  Eqs.  (IV.  10  and 
14),  two  boundary  conditions  are  necessary. 

The  first  condition  comes  from  the  need  to  have  pressure  balance 
across  the  interface.  Integrating  the  radial  component  of  the  equation  of 
motion  across  the  boundary,  we  have 


S'rr 


2^ 


(IV.  15) 


Since  we  require  the  linearized  form  of  Eq.  IV.  15,  we  note  that  the  perturbed 
quantities  are  to  be  evaluated  at  the  unperturbed  boundary  and  the  equilibrium 
quantities,  at  the  perturbed  boundary,  which  has  been  displaced  a distance 
£ = V/y.  We  then  find  the  following  condition  for  pressure  balance 

V 


r = <x 


(IV.  16) 


where  we  have  used  the  equilibrium  pressure  balance  condition  and  the  fact 
that  equilibrium  pressure  was  assumed  uniform. 

The  second  boundary  condition  results  from  the  following  observation. 
Since  the  plasma  is  assumed  to  have  infinite  conductivity  the  electric  field 
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moving  with  the  plasma  must  be  identically  zero,  or 


P V P 

E + - x B =0 
- — -o 

c 

If  we  now  match  fields  in  this  frame;  then,  because  the  tangential  component 
of  electric  field  must  be  continuous  in  any  reference  frame,  we  must  have 

EV  + (V  x B V)  = 0 
-t  o't 

where  the  subscript  t denotes  the  tangential  componert.  This  is  the  second 
boundary  condition.  However  it  can  be  put  in  a more  convenient  form  by 
adding  the  normal  components  of  each  term  to  the  equation  and  taking  the 
curl.  With  the  aid  of  Maxwell's  equations  and  a little  vector  algebra  we 
obtain 


n • 


3B 


n*7  x (V  x B ) 
- - - —o' 


3B 


where  n is  the  normal  to  the  plasma  surface.  Since  — and  V are  zero  in 
equilibrium,  we  can  replace  n by  the  unperturbed  normal  i . Using 
V = we  obtain 


fC  - $ S 


r^oL 


(IV.  17) 


The  final  form  of  the  second  boundary  condition  now  becomes 


(IV.  18) 


We  are  now  in  a position  to  obtain  the  dispersion  relation  for  the  problem  . 
Equations  (IV.  10)  and  (IV.  13)  give  the  result 


(IV.  19) 
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where  prime  denotes  derivative  with  respect  to  argument.  Calculating  the 
perturbed  field  in  the  vacuum  from  Eq.  (IV.  4)  and  substituting  Eq . (IV.  10 
and  19)  into  Eq.  (IV.  16),  we  find 


-A  l'm 


3^ 


r=o l 


, >•  ?>) 


V-7T 


Kw  jko)  C 


h 


O 


(IV.  20) 


where  k*B^  = ~ ^ + kB^  . Then  substituting  from  Eq.  (IV.  19)  and  the 
v 

expression  for  B (from  Eq.  (IV.  14))  into  Eq.  (IV.  18)  results 


i Jk-  _S.V  «0  ijkj) 

p.  (vv  y c^) 


a+  Jt*‘(y-)c=<p 


(IV.  21) 


The  condition  for  non  trivial  solution  of  Eqs.  (IV.  20  and  21)  gives  the 
required  dispersion  relation 

(jk  • iio*-)  i<~ikJ 

#>a  Ik*)  Kd 

+ Jk  (b0vc'»))'  ill**) 


VtT  A a. 


(IV.  22) 


Since  Im' (ka) /im(ka)  >0  and  Km(ka)/Km(ka)<0  it  is  only  the  third  term 
on  the  right  hand  side  of  Eq.  (IV.  22)  which  is  destabilizing.  In  order  to 
see  when  this  term  is  larger  than  the  other  stabilizing  terms,  consider  long 
wavelength  ka<<l.  Using  the  fact  that 


and 


I (ka)/I  (ka)  s-  m/ka 
m m 


K.  (ka) /K  (ka)  w - m/ka, 
m m 
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Eq.  (IV.  22)  can  be  written 


Z 


'/ira'fio  •(IV•23, 


The  most  unstable  k value,  is 


221 

Ck 


■? 


(IV.  24) 


Substituting  this  into  Eq.  (IV.  23),  we  find  that  the  maximum  growth  rate 
is 

(IV.  25) 

so  that  the  m = 1 surface  mode  is  unstable,  m = 0 and  2 are  marginally 
stable  and  all  higher  m (long  wavelength)  are  stable. 

It  is  also  of  interest  to  write  Eq.  (IV.  23)  in  the  form 


A-  j r 


4-rr  fo 

_ _22L  (T£*COk)) 

: v 


(IV.  26) 


a 

Since  only  the  m = 1 mode  is  unstable  we  can  immediately  obtain  the  range 
of  unstable  wavelengths 

0 < I k I < — °9 


a B 


oz 


For  a torus,  the  minimum  value  of  k is  so  a toroidal  plasma  is  stable  if 

a > 1 

(IV.  27) 


?h»)  s 
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Equation  (IV.  27)  is  the  celebrated  Kruskal-Shaf ranov  limit  on  the  q(a) 
value  of  a toroidal  discharge. 

Let  us  see  what  this  condition  implies  for  a tokamak  plasma.  The 

plasma  and  metal  wall  are  shown  in  Fig.  (IV.  1).  The  value  of  q at  the 

plasma  edge  must  be  greater  than  unity.  In  the  vacuum  region  between 
v 1 

r = a and  r = b,  B „(r)  decreases  as  <■  so  that  in  the  vacuum  region 
oo 

q(r)  increases  as  r2,  so  that  q(b)>(/^)2.  Thus,  depending  on  where  the 

plasma  boundary  is,  the  value  of  q at  the  limiter  must  be  somewhat  greater 

than  unity  in  order  for  the  plasma  to  be  stable  to  gross  kink  modes.  This 

lower  limit  on  q(a)  implies  an  upper  limit  on  current  for  stable  configuration. 
0 I 

Using  BQv(r=b)  = where  I is  in  amps,  B in  Gauss  and  r in  cm,  the  maximum 
current  is 

(IV.  28) 


We  now  turn  to  a discussion  of  the  other  type  of  free  surface  mode, 

that  where  the  plasma  has  a uniform  current  up  to  r = a and  vacuum  outside. 

Also  we  assume  B >>B  , R>>a,  and  treat  the  case  where  B is  constant, 
oz  o oz 

In  equilibrium 

_ cJ  I -p  T? 

O-  ~ c zo  Oae>  (iv.  29) 

so  the  pressure  profile  is  parabolic  aid  vanishes  at  r = a,  since  BQ  'v  r 

0O 

and  J is  constant.  We  will  see  that  this  plasma  is  unstable  for  all  m 
zo 

instead  of  just  m = 1 as  was  the  case  for  a surface  current. 

Before  calculating  the  growth  rate,  let  us  show  physically  why  this 
is  so.  As  discussed  in  the  last  chapter,  a perturbation  is  most  likely 
to  be  unstable  if  k»BQ  = 0 since  in  this  case,  it  does  not  couple  to  stable 
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shear  Alfven  waves.  Right  at  the  surface  of  the  plasma,  k*B  = kB  + 

- -o  oz 

mBV  (a) 

od K , and  as  we  just  saw,  the  plasma  is  unstable  when  this  nearly 
r 

vanishes.  However  if  k*BQ  is  zero  at  the  plasma  edge,  it  is  equal  to 

kBzo  inside  the  plasma.  Thus,  modes  which  are  flute  like  (that  is 

perpendicular  to  B ) on  the  plasma  surface  are  not  flute  like  in  the 

plasma  interior,  and  there  is  a strong  coupling  to  shear  Alfven  waves  there. 

This  is  the  explanation  of  the  stabilizing  first  term  on  the  right  hand 

side  of  Eq.  (TV.  22).  Clearly  the  thing  responsible  for  this  stabilizing 

effect  is  the  fact  that  the  magnetic  field  abruptly  changes  direction  as 

one  crosses  the  plasma  boundary.  However,  if  the  axial  current  density 

is  uniform,  the  poloidal  magnetic  field  is  a continuous  function  of  radius 

so  that  a perturbation  which  is  flute  like  on  the  plasma  surface  will  be 

flute  like  in  the  plasma  interior  also.  Thus  we  expect  the  plasma  with 

uniform  current  to  be  less  stable  than  the  plasma  with  a surface  current. 

We  now  proceed  to  derive  the  properties  of  unstable  modes  in  a plasma 

with  a uniform  axial  current  with  B >>B.  and  ^r<<1,  the  so  called 

tokamak  ordering.  Since  the  axial  field  is  very  large,  one  expects  that 

the  most  unstable  motion  could  not  compress  this  field,  implying  two 

dimensional  motion  in  the  r0  plane,  i.e.,  v = B =0.  We  will  assume 

z z 

this  to  be  true  for  now  and  prove  it  a postiori.  The  perturbed  magnetic 

field  in  the  vacuum  is  unchanged  and  is  given  by  Eq.  (IV.  14). 

The  problem  now  is  to  derive  an  equation  for  the  magnetic  field  inside 

the  plasma.  To  do  so,  take  the  z component  of  the  curl  of  the  equation 

of  motion.  Using  the  fact  the  Jq  is  constant  in  the  plasma,  V«J  = V*B  = 0 

and  B is  constant,  we  find 
oz 
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= 


(IV.  30) 


_ 2- 


Relating  J to  B via  Maxwells  current  equation,  and  using  the  fact  that 
z z 

O R.  J-2-r*r+!?Va-V.Vf£r+lr*i?  + MO 

V * P - r r (xv.  : 


we  find 


Ai 
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‘ F / X A ^(rBr)  - i^5r 

yST(  r^r 


(IV.  32) 


where 


B, 


0£> 
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•n  U2« 


-K 


5 M ■ B, 


(IV.  33) 


Note  that  for  the  case  of  constant  current  density  Bqq  is  proportional  to 

r so  F is  constant.  The  other  relation  between  B and  V comes  from 

r r 

Eq.  (IV.  11).  Taking  the  radial  component,  we  find 

y 6>- 


Vr-  = 


(IV.  34) 


Inserting  Eq.  (IV.  34)  into  Eq.  (IV.  32),  we  find  the  single  equation  for 
B 


(IV.  35) 


where  we  have  assumped  p constant. 
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The  factor  outside  is  simply  the  dispersion  relation  for  shear  Alfven  waves 
and  cannot  vanish  for  positive  y.  The  quantity  in  the  square  brackets  must 
then  be  zero.  Trying  a solution  = r11  and  insisting  that  be  well  be- 
haved at  the  origin,  we  find 


A r 


■xn  - 1 


(IV.  36) 


Of  course  A above  and  C from  Eq.  (IV.  14)  can  be  simply  related  by  requiring  that 
B^  be  continuous  across  the  free  surface.  (Since  the  unperturbed  fields  are 
continuous  across  the  free  surface,  all  we  need  do  is  say  that  the  perturbed 
radial  field  is  continuous  across  the  unperturbed  free  surface.)  Using  Eqs. 
(IV.  31  and  34)  we  also  find 


Vs=  L Vr 

V- 


The  dispersion  relation  then  follows  from  assuming  pressure  balance  across 
the  free  surface,  or 


'■2.0 


&.BJ  . 

^90  t )e  ^ -p 

hr 

(IV.  38) 


where  we  have  made  use  of  the  fact  that  B^  = 0 and  that  { Bq  j 


is  continuous. 


In  the  vacuum,  one  can  immediately  derive 


t_ 

r 


Jk  yJ  ( •*») 


(M 

b) 


(IV.  39> 
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It  remains  only  to  calculate  p in  terms  of  B^.  This  comes  from  the  8 
component  of  the  momentum  equation 

? * ±['Jr  8»z+  B.J, 


(IV.  40) 


Now 


Jr-  % v>e 


A 


Sr 


(IV.  41) 


where  we  have  used  Eqs.  (IV.  37b  and  36).  Thus  by  using  Eq.  (IV.  41)  and 
the  expression  for  V~  in  terms  of  from  Eqs.  (IV.  37a  and  34),  Eq.  (IV.  40) 
relates  p to  Br-  Plugging  into  Eq.  (IV.  38),  we  find  the  dispersion  relation 

- (i  jk.  -El-  ( i - V) 

y = ^ ^ Jltpo  \ )r  kn,  (Jea.)  J ^ (!V . 42) 


Since  K^/K^  0,  the  term  proportional  to  F is  stabilizing.  The  term 

7rT  v,  ”^>3* 

proportional  to  F can  be  destabilizing  if  F >0,  or  — p-  ^ 


Clearly  the  plasma  can  be  unstable  in  the  limit  of  F -*  0.  Since  Bqz  is 
constant  and  B^  decreases  as  a function  of  r in  the  vacuum  region,  the 
plasma  can  be  unstable  if  the  singular  surface,  F = 0,  falls  just  outside 
the  plasma.  If  the  singular  surface  falls  inside  the  plasma  F(r  = a)  < 0 
so  the  plasma  is  stable.  Also  if  the  singular  surface  falls  too  far  outside 
the  plasma,  the  second  term  on  the  right  of  in  Eq.  (IV.  42)  dominates  the 
first  and  the  plasma  again  is  stable. 

Let  us  see  what  this  means  for  the  current  buildup  of  a tokamak  plasma. 
As  the  current  rises  the  value  of  q(a)  decreases.  Thus 


F- 
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starts  out  negative  but  as  the  current  builds  up  it  decreases.  When  the 
current  reaches  the  value  at  which 


77 


(IV.  43) 


F changes  from  negative  to  positive  and  the  plasma  is  unstable  to  modes  with 
values  of  m/n  given  in  Eq.  (IV.  43).  However  as  the  current  increases,  the 
mode  becomes  stable  because  of  the  stabilizing  effect  of  the  second  term 
on  the  right  of  Eq.  (IV.  42).  Thus  only  certain  values  of  total  current 
give  instability.  In  building  up  tokamak  current,  the  idea  then  is  to 
program  its  rise  so  that  it  passes  through  the  unstable  values  very  quickly. 
However  the  maximum  plasma  current  will  always  still  be  limited  by  Eq.  (IV.  28). 

We  close  by  determining  the  conditions  under  which  the  motion  is  two 
dimensional,  as  we  have  assumed.  The  velocity  in  the  z direction  comes  from  the 
A component  of  the  momentum  equation.  Assuming  =0,  we  find 

* p.  Vs  --  t-fe  f 

Inserting  for  p from  Eq.  (IV.  40)  and  taking  the  maximum  value  of  y > from 
Eq.  (IV.  42),  we  find 


so  the  motion  basically  is  two  dimensional  as  long  as  kr/2  1. 
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V.  Gravitational  (g)  Modes  in  Slab  Geometry 

This  chapter  discusses  modes  driven  by  a gravitational  force  in 
slab  geometry.  If  gravity  is  present,  the  MHD  equations  for  a perfectly 
conducting,  incompressible  fluid  are 


¥ 

+ 

y v = - v-p  + 

-n  ( 

U ' 

VxVx  & 

Cc) 

V*5 

i?  r 

V-v*  o 
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All  variation  is  assumed  to  be  in  the  x direction. 

At  this  point  we  will  digress  briefly  to  consider  what  the 
significance  of  g is.  First,  and  most  obvious,  g could  correspond  to 
a real  physical  force.  Consider  for  instance  a 9 pinch  implosion  shown 
in  Fig.  VI.  Clearly,  as  the  theata  pinch  implodes,  the  plasma  accelerates 
inward.  In  analyzing  the  local  MHD  stability  of  the  imploding  pinch, 
one  could  work  in  the  reference  frame  in  which  the  fluid  is  locally  at 
rest.  The  inward  acceleration  is  then  described  by  an  outward  inertial 
force.  Also  if  a 0 pinch  plasma  is  rotating  with  angular  velocity  JJ, 
one  could  analyze  its  stability  by  working  in  a rotating  reference  frame. 
Then  the  inertial  force  is  the  sum  of  the  centrifugal,  ft2rir  and  coriolis, 
2ftxV.  The  former  is  outward  and  acts  like  a gravitational  force. 
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Less  obvious,  certainly  less  precise,  but  probably  more  useful, 
the  acceleration  g can  be  used  to  model  complicated  geometric  aspects  of 
the  magnetic  field.  To  see  this,  let  us  examine  the  right  hand  side  of 
Eq.  (lb)  without  a g term 

-vP+ £&■?)>§  *- *?- «*-  (*■•)  ! 

J 

In  slab  geometry,  the  lines  of  B in  equilibrium  are  assumed  to  be  straight. 

Therefore  any  equilibrium  variation  in  B"  is  balanced  by  a change  in 
pressure.  (That  is  (B-V)B  = 0 in  slab  geometry  in  equilibrium).  In  a 
more  realistic  equilibrium,  there  are  additional  forces  on  the  fluid 
which  result  from  geometric  aspects  of  B.  For  instance  in  cylindrical 
geometry,  the  lines  of  force  are  curved.  In  toroidal  geometry,  not  only 
are  the  lines  of  force  curved,  there  is  also  a change  in  the  magnitude 
of  the  toroidal  field  resulting  from  the  (geometric  effect  of  the)  decrease 
of  field  with  major  radius.  The  question  is  whether  these  additional 
forces  on  the  fluid,  which  arise  from  geometric  dependence  of  B,  can  be 
modeled  with  a gravitational  force. 

Of  course  there  is  clearly  no  way  in  \hi.ch  this  simple  gravitational 
acceleration  can  model  all  of  the  very  complicated  effects  of  field 
inhomogeniety  and  curvature.  However,  as  we  will  see,  in  a collisionless 
plasma  at  least  (collision  frequency  much  less  than  ion  cyclotron  frequency) 
one  can  convincingly  argue  that  many  of  the  important  effects  of  complicated 
field  structure  can  indeed  be  approximately  be  modeled  by  gravity. 

The  force  on  a fluid  element  of  course  is  just  the  sum  of  forces 
on  all  individual  particles  comprising  the  fluid  element.  If  the  field 
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lines  are  curved,  having  radius  R,  the  particles  feel  an  outward 

centrifugal  force  MV2/R  where  V is  the  velocity  parallel  to  the  line  of 

il  J 

force,  and  R-1  = -(B*V)B/B2.  Thus  the  force  on  the  particle  resulting 
from  the  curvature  of  the  field  line  is 


- - 


Mv‘  fl-  V B 


In  addition  to  the  force  resulting  from  the  curvature  of  the 
field  lines,  there  is  also  a force  arising  from  the  gradient  in  magnitude 
of  B perpendicular  to  B.  Consider  for  instance  a magnetic  field  in  the  z 
direction  which  increases  in  the  x direction.  The  force  in  the  negative 
x direction  on  the  right  hand  side  of  the  particle  orbit  is  larger  than 
the  force  in  the  positive  x direction  on  the  left  hand  side  of  the  particle 
orbit,  as  shown  in  Fig.  V2.  Thus  there  is  a net  force  in  the  negative  x 
direction 


U-  - ywi 
•?!?/  ’ 
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(These  average  forces  are  calculated  in  much  more  detail  in  any  standard 
plasma  physics  textbook.  See  for  instance  Schmidt,  Physics  of  Hot  Plasmas, 
or  Krall  and  Trivelpiece,  Plasma  Physics.)  The  equivalent  gravitational 
force  g can  now  be  calculated  from  Eq.  (V  3 and  4).  Of  course  in  using 
Eq.  (V4)  for  Fi,,  one  must  only  use  that  part  of  V j B ] which  is  due  to 
geometric  effects  and  not  use  that  part  which  is  balanced  by  a pressure 
gradient.  This  latter  part,  of  course,  is  perfectly  well  described  in  a 
slab  model. 

For  our  purposes  now,  we  simply  model 

r Vt>o? 
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where  iR  is  a unit  vector  pointing  outward  along  the  radius  of  curvature 

of  the  field  line,  and  is  a typical  thermal  velocity  (that  is  is 

1/2 

roughly  either  the  ion  thermal  velocity  or  else  (Me/ML)  times  the 
electron  thermal  velocity) . 

After  this  digression  on  the  significance  of  gravity,  we  return 
to  a discussion  of  g modes  as  described  by  the  linearized  version  of 
Eq.  (VI).  All  variation  is  taken  to  be  in  the  x direction  and  g also 
points  in  the  x direction.  The  magnetic  field  points  in  the  z direction 
and  is  taken  to  be  uniform.  The  basic  g mode  is  then  very  easily  derived. 
Assuming  a time  dependence  exp  yt  for  perturbed  quantities  the  linearized 


versions  of  Eq.  (1)  are 


+ Vv  *© 
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where  now  perturbed  quantities  have  no  subscript  and  unperturbed  quantities 
have  a subscript  zero.  Let  us  further  simplify  by  assuming  a spatial 
dependence  of  exp  i ky  (that  is  no  x or  z dependence)  for  perturbed 
quantities.  Then  inserting  from  Eq.  (V6a)  into  Eq.  (V6b)  and  taking  the 
x component,  we  find  ^ 

If1--  - (%?) 

r- 

Thus,  if  g and  3po have  opposite  signs  the  plasma  is  unstable.  These 
are  the  classic  instabilities  one  has  if  magnetic  field  lines  are  convex 
to  the  plasma  like  the  classic  magnetic  mirror  arrangement  shown  in 
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Fig.(V3).  Here,  the  quantity  g as  approximated  ly  Eq.  (V5)  is  outward,  or 


in  the  positive  radial  direction.  The  density  gradient  however  is  in  the 


negative  radial  directs  1 so  that  g 


< 0 and  the  plasma  is  unstable. 


Notice  also  that  this  instability  is  extremely  violent.  Assuming  that 
L = (p  '"o  | is  of  order  of  R which  is  of  order  of  the  size  of  the 


L” 5 (°° 


plasma  R,  the  growth  time  is  of  order  R/V^.  That  is  the  time  for  the  mode 
to  grow  is  roughly  the  time  for  the  plasma  to  escape  from  the  mirror  trap 
traveling  at  the  thermal  velocity;  in  other  words,  there  is  no  confinement 
of  the  plasma  at  all! 

Clearly,  in  order  to  confine  the  plasma,  some  means  must  be  devised 

to  stabilize  these  modes.  One  possibility  is  to  utilize  a cusp  (like  minimum  B)  type 

configuration  like  that  shown  in  Fig.  (V4)  so  that  g > 0 everywhere. 

3x 

Another  possibility  is  to  find  some  way  to  stabilize  these  modes.  One  way 

to  stabilize  them  is  by  magnetic  shear.  A discussion  of  this  stabilization 

mechanism  will  occupy  the  remainder  of  this  chapter.  In  order  to  set  the 

stage  for  our  calculation  of  shear  stabilization  of  g modes,  we  begin  by 

considering  these  modes  when  k^  r 0. 

If  perturbed  quantities  vary  as  exp  i (k  y + k Z) , tl  n the  x 

y z 


component  of  Eq.(V6b)  becomes 

_ L >2  B, 
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To  obtain  B..  in  terms  of  V , one  uses  the  x component  of  Eq.  (Vic). 
Inserting  for  B^  in  terms  of  V , we  find  that  the  dispersion  relation  is 


y = 


- >z  Bfl 
V-rr  po 
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(X  1) 


51 


3p0 

Notice  that  for  k 2V  2 > g -r — /p  the  mode  is  stable.  Clearly  as  k 

2 A dX  0 Z 

increases,  the  unstable  g mode  attempts  to  couple  to  a stable  shear  Alfven 
mode.  Postulating  a k^  sufficiently  large  that  y2  given  by  Eq.  (V9)  is 
negative  is,  of  course,  not  equivalent  to  stabilizing  the  system.  In  an 
unsheared  field,  the  plasma  is  still  free  to  pick  a parallel  wave  number 
equal  to  zero.  The  basic  idea  behind  shear  stabilization  is  to  remove 
this  freedom  by  imposing  an  equilibrium  variation  on  the  direction  of  B. 

Let  us  now  say  that 

B.m  = 73*  (i*  * £z)  (y: 10 ) 

where  we  imagine  x <<  L . Since  this  additional  x variation  is  imposed 
on  the  system,  all  perturbed  quantities  have  the  functional  form 

f (?) 

where  we  have  set  k =0.  Notice  however  that  the  variation  in  a 

z 

perturbed  quantity  is  no  longer  perpendicular  to  B . Indeed 

fa-i)  *}***)*  ‘J  s’£  &n) 

Thus  the  presence  of  shear  forces  any  perturbed  quantity  to  vary  parallel 
to  Bq  as  long  as  x a 0.  This  forced  variation  parallel  to  B will  tend  to 
couple  to  shear  Alfven  modes  and  stabilize  the  plasma.  We  now  proceed  to 
calculate  just  how  much  shear  is  needed  for  stabilization. 

In  a sheared  field,  the  equation  for  the  perturbed  x velocity  is 
slightly,  but  not  much  more  difficult  to  specify.  Because  of  the  inherent 
x dependence  of  Bq,  one  can  no  longer  assume  that  perturbations  are 
independent  of  x.  To  begin,  take  the  curl  of  the  perturbed  momentum  equation 
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where  we  have  made  use  of  the  fact  that  B=—  VxVxB  . Now  consider 

~ y - - -o 

the  x component  of  Eq.  (V  12).  Recalling  that  = 0,  a short  calculation 
gives  the  result 

[/<■>'■  (I,3> 
Equation  (V  13)  above  dictates  that  either  = 0,  or  else  the  quantity  in 
the  square  brackets  vanishes.  The  latter  case  implies  shear  (stable) 

Alfven  waves  at  each  point  x.  Since  these  are  not  unstable,  we  assume 
the  former,  ie 

V = 0 (V  14) 

z 

From  the  fact  that  = 0,  one  can  show  easily  from  Eq.(V  Id)  that  also 

B = 0 (V  15) 

z 

Thus  the  perturbation  is  two  dimensional  in  the  x-y  plane.  The  incompressibility 

condition  then  relates  V to  V and  B to  B 

X y x y 

Then,  taking  the  z component  of  Eq.  (V  12)  and  making  use  of  Eq.  (V  16) 
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above,  we  find  the  following  simple  equation  for  V : 


Notice  that  the  term  in  the  parentheses  on  the  right  hand  side  is  simply 
the  local  dispersion  relation,  Eq.  (V  9),  if  one  assumes  k2  = k2(x/Lg)2. 
However,  this  is  no  longer  the  dispersion  relation  in  the  sheared  system. 
For  instance  if  Eq.  (V  9)  were  initially  satisfied,  disturbances  near 
x = 0 would  grow  while  those  far  from  x = 0 would  not.  A strong  x 
dependence  would  then  be  induced  by  the  local  growth,  and  eventually  the 
first  term  of  Eq.  (V  17)  would  also  be  important.  The  problem  then  is  to 
find  the  eigenfunctions  and  eigenvalues  of  Eq.  (V  17).  In  the  limit  of 
x °°,  the  equation  approximately  reduced  to 


with  solution 


&I7  ) 


clearly  the  proper  boundary  condition  is  that  approaches  zero  at  both 


x = ± a°.  Since 


P‘ 


y(rj£  = <?>* 


*/rr 


as  long  as  y2  > 0,  Eq.  (V  17)  has  no  singular  points  on  the  real  x axis. 
Therefore  one  could  easily  solve  this  equation  numerically  for  the 
eigenvalues  y2. 

We  now  see  what  insights  may  be  obtained  analytically.  One  can 
immediately  show  that  if  g > 0,  no  solution  with  y2  > 0 can  satisfy 
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the  boundary  conditions.  To  do  so,  note  that  py2  + k^y25-)2  - — £ Q(x)  is 

L 4tt 
s 

then  greater  than  zero.  Assume  that  at  x 4 - °°,  V is  well  behaved  and 

x 

also  that  V >0.  Then 

where 

G(x)=  Q(*-)+3y?  &•?*) 

3p  3V 

0 X 

which  is  greater  than  zero  since  g . Therefore  - — > 0 for  all  x 

3x  3x 

so  that  the  solution  for  is  a monotonically  increasing  function  of  x. 
Hence,  it  does  not  satisfy  the  proper  boundary  condition  as  x -*•  + <=°. 

Clearly,  in  order  for  Vx  to  satisfy  the  proper  boundary  condition  as 
x ->•  + <*>,,  there  must  be  some  region  of  x where  G(x)  < 0. 

Now  let  us  examine  how  an  unstable  eigenfunction  can  be  constructed. 

3P0 

If  g -tT  is  greater  than  zero  the  function  is  monotonically  increasing, 

and  as  x it  diverges  as  exp  kx.  Now  depending  upon  the  sign  of  G(x) 

(Eq.  V 22),  the  equation  has  two  fundamentally  different  types  of  behavior. 

If  G(x)  > 0,  the  solution  for  is  a monotonically  increasing  function 

of  x as  we  have  just  discussed.  On  the  other  hand,  if  G(x)  < 0,  the 

solution  for  is  generally  oscilliatory  in  x.  The  local  wave  number 

1/2 

is  given  very  roughly  by  k(-G(x)/Q(x))  . If  G(x)  changes  sign  as  x 

varies,  we  will  denote  the  regions  of  G(x)  > 0 as  monotonic  regions  and 

the  regions  where  G(x)  < 0 as  oscilliatory  regions.  We  have  just  shown 

that  the  regions  x -*■  ± 00  are  monotonic. 

If  g 3po  < 0,  there  is  an  oscilliatory  region  specified  by 
^3x 
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Since  the  wave  number  of  the  oscillation  in  the  oscilliatory  region  is 

1/2 

given  roughly  by  k(-G(x) /Q(x) ) , it  is  clear  that  increasing  y increases 

the  wavelength  in  .he  oscilliatory  region.  Let  us  imagine  that  a y = y_ 

can  be  found  such  that  V = 0 at  x = x . Then,  since  x > x is  a 

x os  os 

monatonic  region,  the  solution  for  for  x > x^  diverges  toward  minus 

infinity  as  x -*■  + °°  as  shown  in  curve  A in  Fig.  (V  5) . As  y is  increased, 

the  wavelength  in  the  oscilliatory  region  increases  so  that  the  solution 

for  Vx  will  not  curve  downward  as  much  here.  Thus  for  a sufficiently 

large  y 5 y+,  the  solution  for  diverges  toward  positive  infinity  as 

x+»>as  shown  in  Fig.  (V  5)  curve  B.  Clearly,  there  exists  some  y - YQ 

where  y c yQ  < y+  for  which  the  solution  to  Eq.  (V  17)  assymptotes  to 

zero  for  through  positive  values  of  V+  as  shown  in  Fig.  (V  5)  curve  C. 

This  then  is  an  unstable  eigenfunction  with  eigenvalue  y . 

One  can  now  deduce  the  following  criterion  for  determining  whether 

Eq.  (V  17)  has  unstable  eigenvalues.  First  select  a y and  assume  Vx  is 

well  behaved  as  x - 00  and  integrate  the  equation  towards  plus  infinity. 

If  there  is  one  zero  of  V before  the  solution  diverges  at  x = + <*>,  then 

x 

there  will  be  one  eigenfunction  with  a growth  rate  larger  than  this  assumed 
y.  Similarly,  if  there  are  two  zeros  of  before  the  solution  diverges 
for  x -*■  + °°,  there  are  two  eigenfunctions  with  growth  rate  larger  than 

the  assumed  y;  and  so  on.  Clearly,  the  more  nodes  to  the  eigenfunction, 
the  smaller  is  y.  The  equation,  and  this  technique  of  determining  a sequence 
of  eigenvalues  is  closely  related  to  the  Sturm  - Liousville  equation 
described  in 
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Morse  and  Feshback,  Methods  of  Theoretical  Physics,  Vol.  I,  Chapter  VI 
or  almost  any  advanced  calculus  textbook. 


Now  it  is  possible  to  derive  very  quickly  an  approximate  criterion 

for  shear  stabilization  of  g modes.  In  order  that  V (x  ) % 0,  it  is 

x os 

necessary  that  there  be  roughly  one  half  a wavelength  of  oscillation  in 
the  oscilliatory  region.  Approximating  the  wave  number  of  this  oscillation 
by  the  square  root  maximum  value  of  k2G  divided  by  the  square  root  of  the 
maximum  value  of  Q,  we  find  instability  if 


For  Lg  much  larger  than  specified  by  Eq.  V 24,  the  mode  grows  with  growth 
rate  given  roughly  be  Eq.  (V  7). 

Equation  (V  25)  is  then  a condition  for  shear  stabilization  of  g modes 
which  have  scale  length  of  order  x . One  might  think  that  by  localizing 

OS 

the  mode  very  close  to  x * 0,  and  by  greatly  decreasing  the  growth  rate, 
the  wavelength  of  the  oscillation  for  V in  the  oscilliatory  region  could 
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be  significantly  reduced  so  that  unstable  eigenfunctions  could  be 
constructed  even  if  Eq.  (V  25)  is  violated.  This  possibility  will  now 
be  examined. 


Y 
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For  very  small  x and  y,  the  equation  for  can  be  written  as 


2_ 
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In  Eq.  (V  26)  above,  the  y2  is  retained  in  the  left  hand  term  so  that  the 


solution  is  non  singular  at  x = 0.  However  clearly,  for 

z 


2 » L 


1 V-r 
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the  y term  in  Eq.  (V  26)  is  not  important.  Therefore,  we  can  neglect  y, 
but  restrict  ourselves  to  the  range  of  x denoted  by  Eq,  (V  27).  By  consid- 
ering the  limit  x + 0 we  can  see  how  the  solution  behaves  for  small  x and 
y,  and  particularly,  we  can  see  whether  the  solution  has  a node.  (Naturally 
one  cannot  connect  the  solution  at  positive  x to  the  solution  at  negative  x 
in  this  way,  because  to  do  so  one  must  also  have  the  solution  in  the  region 
where  Eq.  (V  27)  is  violated.) 

If  y =0  in  Eq.  (V  25),  one  can  solve  by  assuming  Vx  = x11.  (Appendix 
is  a review  of  some  aspects  of  power  series  solutions  of  second  order 
differential  equations).  In  this  case  one  can  solve  for  n and  find 

i (>+  **)*  &"*') 

where 

*3  Jr  £ /tl  fa:  ) 

The  condition  for  the  exponent  n to  be  complex  is  simply 

9 < -Jf 


58 


(T3c) 


However  if  n is  complex,  n = nf  + in^  the  behavior  of  Vx  is  oscilliatory, 
and  the  wavelength  of  oscillation  approaches  zero  as  x -*■  0,  that  is 

Vr  - z ^ x) 

Thus  there  are  an  infinity  of  nodes  of  V in  the  limit  of  x 0.  This 
means  that  there  are  an  infinity  of  unstable  eigenfunctions.  The  growth 
rates  of  these  eigenfunctions  approach  zero  as  the  number  of  nodes  increase. 
The  condition  for  such  localized  unstable  modes  is  given  by  Eq.(V  30), 
while  the  approximate  condition  for  stabilizing  modes  spread  over  x is 

given  by  Eq.  (V  24).  Notice  that  the  equations  have  roughly  the  same 
form  except  that  about  three  times  as  much  shear  (that  is  Lg  smaller  by 
one  third)  is  required  to  stabilize  modes  localized  over  x^  than  is 
required  to  stabilize  modes  localized  right  at  the  position  x = 0.  A 
condition  for  the  occurrence  of  unstable  modes  is  then  given  by  Eq.  (V  30). 
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tiiavorable  curvature 


Fig.  V 5 — Diagram  of  how  to  form  the  eigenfunction  and  calculate  the  eigenvalue  (7)  of  Eq.  (V  17). 
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VI . Resistive  g Modes 

In  the  previous  chapter  it  was  shown  how  shear  in  the  magnetic  field 
produces  a stabilizing  effect  on  instabilities  driven  by  a gravitational 
field  g.  In  this  chapter  we  shall  again  consider  the  stability  of  a plasma 
L supported  against  gravity  by  a sheared  magnetic  field,  but  with  the  addition 

of  a further  physical  effect,  namely  resistivity.  For  the  g modes  of  the 
previous  chapter  the  ideal  MHD  model  was  used  for  which  the  field  and  fluid 
are  frozen  together.  However,  for  the  instabilities  to  be  described  in 
this  chapter,  the  magnetic  field  and  fluid  become  decoupled  due  to  the 
presence  of  resistivity.  It  will  be  found  that  the  instability  is  again 
localized  at  the  resonant  surface  where  k*B  = 0,  which,  as  we  will  see,  is 
the  only  place  where  a small  resistivity  can  introduce  a significant  amount 
of  decoupling  in  motion  of  field  and  fluid.  The  net  effect  of  resistivity 
on  g modes  is  to  nullify  the  stabilizing  effect  of  magnetic  shear.  An  ideal 
g mode  which  was  stabilized  by  the  presence  of  magnetic  shear  will  be  shown 
to  be  unstable  - albiet  at  a reduced  growth  rate  - when  the  presence  of  a 
small  amount  of  resistivity  is  allowed  for.  Thus,  if  the  effect  of  magnetic 
shear  is  compared  to  a dyke,  then  resistivity  turns  the  dyke  into  a leaky 
one. 

Let  us  now  describe  the  resistive  g mode  in  slab  geometry,  as  in  Chapter  V. 
The  magnetic  field  is  taken  as 

BJ*)  = (z  30  fl 

where  B^  is  a constant.  The  gravitational  acceleration  is  in  the  x direction 
and  all  equilibrium  fluid  quantities  also  vary  in  only  the  x direction.  We 
treat  the  plasma  as  an  incompressible,  weakly  resistive  fluid.  The  linearized 
equations  of  motion  are 
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X7  * /Vx  *6) 
y^r  v 


v-v  = o 


£n.*0 

(^r  i) 
Cozr.  v ) 


Note  that  we  have  treated  the  resistivity  as  uniform  in  obtaining  Eq.  (VI.  2). 
Since  the  equilibrium  varies  only  with  x,  we  know  that  the  eigenfunction  of 
all  perturbed  quantities  varies  as  f(x)  exp  (yt  = iky).  We  then  find  the  V 

z 

and  decouple  from  the  remaining  variable  so  that  only  the  x and  y components 

of  Eq,  (VI.  1 and  2)  are  required.  (The  variables  V^  and  B^  describe  the 
properties  of  shear  Alfven  waves  so  that  we  may  put  both  of  these  variables 
equal  to  zero  when  analyzing  the  behavior  of  linearly  independent  resistive 
g modes) . 


The  x and  y components  of  Eq.  (VI  1 and  2)  give 


>/>.  v,  --  |f  ^ ' sir  &(%£ £zr  ^ 
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We  now  simplify  the  problem  by  expressing  all  variables  in  terms  of  either 


V or  B . Using  Eqs.  (VI.  6 and  4)  we  find 

O - y /°0  V*  _ t Bo  y? 

Substituting  Eqs.  (VI.  3 and  9)  into  Eq.  (VI.  5)  and  using  Maxwell's 
equation  V*B  = 0,  gives 


frZE  7 ) 


2. 


Y-57" 

(yr  )o 


We  must  now  solve  the  pair  of  equations  (VI.  7)  and (VI.  10).  Notice 
that  with  n = 0 in  Eq.  (VI.  7),  the  equations  describing  g modes  in  ideal 
MHD  results. 

Before  solving  these  equations  with  n 4 0,  let  us  consider  the  effect 
of  resistivity  on  the  stability  of  the  system.  As  already  mentioned,  the 
instability  will  be  localized  around  the  resonant  surface  k*BQ  = 0.  This 
is  also  the  position  where  the  motion  of  the  fluid  and  field  can  most  easily 
decouple.  The  form  of  Ohm's  law  which  has  been  used  to  give  Eq.  (VI.  2)  is 

ys=  € ?.  (yr  n) 

For  q = 0 the  condition  of  frozen  in  field  lines  is  simply 

/ = - -jr-fi  ?. 

On  the  other  hand,  the  opposite  extreme  to  this  is  clearly  when  the  fluid 
moves  freely  through  the  field  lines,  i.e.,  V ^ 0 but  E « 0.  In  this  case, 

Eq.  (VI.  11)  gives 


c is 


% y 
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One  part  of  the  force  exerted  on  the  plasma  is  then 

3--,  • f.X\z  Vr 

-==-  / - V D„  / 


F- 


£ST  H ) 


From  the  negative  sign,  we  see  that  the  force  is  a restraining  force  and 
furthermore,  for  x ^ 0 and  small  n,  this  force  is  very  large.  This  only 
means,  of  course,  that  the  plasma  does  not  slip  through  the  field  lines, 
rather  the  field  remains  frozen  into  the  flow.  However  for  x = 0 (that  is 
k-B  = 0),  the  restraining  force  becomes  small  and  the  plasma  can  leak 
through  the  field.  We  therefore  see  the  importance  of  the  resonant  surface 
and  the  reason  that  resistivity  can  be  important  at  this  point,  but  un- 
important everywhere  else  in  the  plasma. 

Let  us  assume  that  in  the  viscinity  of  the  resonant  surface  the  plasma 
flow  is  decoupled  from  the  field,  so  that  Eq.  VI.  7 becomes 


As  we  will  see  shortly,  one  simplifying  feature  of  resistive  g modes,  as 
compared  to  for  instance  tearing  modes  (next  chapter),  is  that  resistive  g 
modes  are  localized  about  the  resonant  surface.  Therefore,  for  resistive  g 
modes,  (Eq.  (VI.  5)  can  be  taken  to  apply  everywhere. 


Substituting  Eq.  (VI.  15)  into  Eq.  (VI.  10),  we  obtain 

1 2 i 


1 P/ 


2* 


2^ 


V* 


tszr/0 


The  problem  is  now  reduced  to  the  solution  of  a single  second  ordinary 
differential  equation,  the  quantum  mechanical  harmonic  oscillator  in  fact. 
Since  we  are  interested  in  solutions  localized  around  x = 0,  we  neglect  the 
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in  Eq.  (VI.  5). 


3ZB 

, 

dependence  of  p on  x and  also  neglected  kzB  compared  to  Jx2 

o x 

Finally  we  have  assumed  that  the  growth  rate  y is  much  less  than  the  growth 

3p0\l/2 

rate  in  an  unsheared  field  ( g/p  ).  Introducing  the  constants 


(g/co  sr)- 


fi,  y 7 A 4 


and 


s--  - 

3 «£•  v ±i_ 

T & / a1 

Eq.  (VI.  16)  can  be  rewritten 

A J'Vx  j. 

* 

(S  - = <= 

Introducing  a change  in  variables 

/ 

X-  > 

A x 

Eq.  (VI.  8)  becomes 

( fa  - 

[SL  17°-) 
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The  solutions  to  Eq.  (VI.  19)  which  vanish  at  x = ± 00  are  the  Hermite  functions 


V 

and  the  eigenvalues  are 


By  = + I 

S' ft'/* 


('&JL 

[yr.  *i) 


Notice  that  there  can  only  be  a real  eigenvalue  y,  for  B > 0.  From  the 
definition  of  B,  Eq.  (VI.  17. b),  this  means  gPQ'<0;  that  is  density  gradient 
is  opposite  to  g so  that  the  g mode  in  an  unsheared  field  is  unstable. 

Equation  (VI.  21)  is  the  dispersion  relation  for  the  resistive  g mode. 
Substituting  for  A and  B,  we  find 


Y- 


Jnt!  ( </  />*  S. 
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This  shows  the  characteristic  dependence  of  the  growth  rate  of  the  resistive 

g mode  on  the  one  third  power  of  resistivity.  Notice  also  that  this  mode 

cannot  be  removed  by  shear.  As  Lg  decreases  (increasing  shear),  the  growth 

rate  is  reduced  but  there  is  no  critical  value  at  which  it  vanishes.  Hence 

the  earlier  analogy  with  a "leaky  dyke". 

Finally,  let  us  use  the  above  solution  to  check  the  validity  of  the 

approximation  which  enabled  us  to  carry  out  this  analysis.  This  was  the 

2 32B 

nc4-  v 

neglect  of  the  yB  compared  with  . v • in  Eq.  (VI.  7).  This  approximation 

x 4 it  dr 

is  only  expected  to  be  valid  over  a small  distance,  (say  xc)  where  is 
given  approximately  from  Eq.  (VI.  16)  by  the  condition 


Using  this  to  compare  the  orders  of  magnitude  of  the  neglected  term  yB 


From  the  expression  for  g given  in  the  previous  chapter 

F.  P. 

a v rJ  __ 

3 A,  (ff  /*  r (31  i-3. 


/>.  r (V:  -B.*; 


where  R,  the  radius  of  the  field  line  is  r(B2  + B2)/F“,  and  L v Rq  (R  is 

z 8 0 s 


the  uiajor  radius  of  the  torus),  the  condition  for  validity  becomes 

■A  f « i f- 

s;. ' s.\ 

This  condition  is  well  satisfied  in  all  current  tokamaks  and  pinches 


(-2E  wj) 
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VII.  The  Tearing  Mode 

In  this  chapter  we  turn  to  a discussion  of  one  of  the  most  important, 

and  interesting  modes  in  all  of  MHD  theory,  the  tearing  mode.  This  mode  is 

quite  different  from  the  resistive  g mode.  There,  resistivity  allowed  the 

plasma  to  slip  through  the  sheared  field,  rather  like  (but  much  slower 

than)  the  case  where  there  was  no  shear.  The  tearing  mode  is  a totally 

different  type  of  fluid  motion  and  it  is  allowed  only  because  of  the 

presence  of  resistivity.  In  the  resistive  g mode,  the  fastest  growing 

mode  had  V an  even  function  of  x.  For  the  tearing  mode  however,  V is 
x x 

an  odd  function  of  x so  that  for  a particular  y,  two  fluid  elements  on 
opposite  sides  of  the  singular  surface  are  flowing  either  towards  each 
other  or  away  from  each  other. 

Before  actually  discussing  the  tearing  mode,  it  is  worthwhile  to 
examine  this  type  of  motion.  By  doing  so,  we  will  show  that  if  magnetic 
field  lines  can  reconnect,  certain  restraining  forces  disappear  so  that 
the  equilibrium  is  more  likely  to  be  unstable.  Let  us  consider  two 
dimensional  incompressible  motion  in  the  xy  plane  where  initially 

b --  % fL  (vm;  I ) 


Then  specify  the  following  velocity  profile  in  the  x direction 


- ya<wi> 
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V*  r 


- V — CasJk 
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where  6 is  taken  to  be  very  small.  Using  the  fact  that  the  flow  is  two 
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dimensional  and  incompressible,  we  find 


O S < * 

O < <^  < X 


( oS) 

( 3nr  3) 


V3(-!0-  VjI*)  (b) 

Thus  in  the  separation  layer  |x|<<$,  there  is  very  rapid  flow  in  the  y 

i 

direction.  The  velocity  stream  lines  are  as  shown  in  Fig.  (VII  1).  The 
two  fluids  on  opposite  sides  of  the  singular  surface  x = 0 stream  towards 
each  other  for  y > 0,  recoil  near  x = 0 and  then  stream  away  from  each 
other  at  some  y < ir/k.  Similarly  two  fluid  elements  streaming  towards 
each  other  at  y < 0 recoil  and  stream  away  somewhere  around  y > - ir/k. 

The  fluid  exactly  at  y = 0 initially,  does  not  stream  away  toward  either 
positive  or  negative  y,  but  piles  into  the  stagnation  point  at  x = y = 0. 

We  will  call  this  type  flow  a tearing  flow  or  tearing  motion. 

Now  imagine  that  the  fluid  is  magnetized  with  the  field  given  by 
Eq.  (VII  1).  In  ideal  MHD,  the  magnetic  field  is  frozen  into  the  flow. 

This  allows  us  to  see  just  how  the  magnetic  field  reacts  to  the  flow 
pattern  shown  in  Fig.  (VII  1).  In  Fig.  (VII  2a)  is  shown  a portion  of 
two  magnetic  field  lines  - ^ < y < ^ initially  at  x and  -x  where 
| x|  >6  at  time  t = 0.  On  one  of  these  field  lines  in  Fig.  (VII  2a)  are 
a series  of  dots  (A-A').  Since  the  field  is  frozen  into  the  flow,  it 
must  always  pass  through  these  dots  no  matter  where  they  go.  Thus  by 
following  the  position  of  these  dots,  we  can  construct  the  field  line 


at  later  times.  The  field  line  for  positive  x can  of  course  be  drawn 
in  an  analogous  way.  Since  V^  = V^  = 0 at  y = ± the  points  A and  A' 
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do  not  move,  while  the  other  points  B-E  move  toward  x = - 6 with  progressively 
higher  velocity.  Thus  at  a later  time  the  field  lines  bulge  as  shown  in 
Fig.  (VII  2b).  Now  consider  a still  later  time  when  points  D,  E and  D'  have 
entered  the  separation  region  - 6 < x < 0.  Point  D shoots  upward,  D'  down- 
ward and  E remains  on  the  line  y = 0,  always  approaching  the  separation 
point  x = y = 0.  Thus  the  field  lines  are  distorted  as  shown  in  Fig.  (VII  2c). 

At  a later  time  this  distortion  is  further  accentuated  and  the  field  line  j 

looks  as  shown  in  Fig.  (VII  2d).  Notice  that  the  field  lines  get  very 
twisted  and  contorted  by  this  type  of  motion. 

This  contortion  of  the  field  line  strongly  affects  the 
fluid  motion.  The  field  line  is  sharply  curved  at  the  points  marked  with 
a * on  the  right  hand  side  of  Fig.s  (VII  c and  d) . However,  as  we  have 
seen  in  Chapter  III,  a curved  field  line  exerts  a tension  along  the  line 
rather  like  a rubber  band.  Thus  the  magnetic  forces  at  the  starred  points 
in  Figs.  (VII  2c  and  d)  (and  of  course  at  the  analogous  points  on  the  left 
hand  graph)  are  trying  hard  to  snap  the  line  back  into  its  original  shape. 

That  is,  the  magnetic  forces  provide  a strong  restoring  force  which  tends 
to  prevent  the  type  of  fluid  motion  given  by  Eq.  (VII  2 and  3). 

Now,  let  us  examine  what  effect  non  zero  resistivity  has  on  this 
motion.  If  the  plasma  has  resistivity,  the  field  lines  can  break  and 
reconnect.  Notice  also  that  in  addition  to  the  sharp  corners  induced 
at  the  stars  in  Fig.  (VII  2c  and  d) , there  are  also  long  lines  of  weak, 
oppositely  directed  fields  forced  right  next  to  each  other  near  the 
singular  surface  x = o.  Thus  if  resistivity,  and  thereby  magnetic 
diffusion  is  allowed,  the  field  lines,  rather  than  looking  like  those  in 
Fig.  (VII  2c  and  d)  would  look  like  that  shown  in  Fig.  (VII  3). 
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Notice  that  the  field  pattern  in  Fig.  (VII  3)  does  not  have  nearly 
as  many  sharp  corners  as  that  in  Fig.  (VII  2 c or  d).  Thus,  it  provides 
a much  weaker  restraining  force  against  the  tearing  motion  described  by 
Eq.  (VII  2 and  3).  Indeed,  as  we  will  see  in  this  chapter,  configurations 
which  are  stable  to  tearing  motion  in  ideal  MHD  can  be  unstable  if  the 
resistivity  is  non  zero. 

We  now  continue  by  writing  out  the  equations  for  the  linear  stability 

of  the  system.  We  assume  that  B = B i+B  (x)iy  where  B = o at  x = 0. 

-oozoy  oy 

The  procedure  is  exactly  the  same  as  in  Chapter  V,  namely  take  the  x 

component  of  the  curl  of  the  momentum  equation  to  show  that  V = 0.  The  z 

z 

component  of  Eq.  (V  lc)  then  shows  that  B^  = 0.  Then  take  the  z component 
of  the  curl  of  the  momentum  equation,  use  the  fact  that  V and  B are  two 
dimensional  and  incompressible,  and  arrive  at  the  result 


i 


(°o 

ik-  ^ * 


+ 3*  I (ye 


assuming,  as  usual  that  perturbed  quantities  vary  as  exp(6t  + iky).  Notice 
that  Eq.  (VII  4)  combined  with  the  x component  of  Eq.  (V  lc)  would  give 
Eq.  (V  17)  except  that  in  the  former  g is  assumed  to  vanish,  while  in  the 
latter,  the  second  derivation  of  the  current)  is  assumed  zero. 

If  the  resistivity  is  present,  the  equation  relating  B to  V is 

* % - ‘ A 8,  v> + g(. 3 - x &nr  ^ 

Let  us  first  prove  that  this  configuration  is  stable  in  ideal  MHD  (i.e.  if 

n = 0) . Substituting  B = ikB  V fy  into  the  right  hand  side  of  Eq.  (VII  4), 

x yo  x n ' * 


Notice  that  for  real  y,  the  right  hand  side  of  Eq.  (VII  6)  is  positive 

for  all  x.  Thus,  as  discussed  in  Chapter  V,  is  a monotonically  increasing 

3 V 

function  of  x for  all  x as  long  as  V and  x are  both  positive  in  the 

-Tx 

limit  of  x -+  °°.  Hence  there  is  no  way  to  have  a solution  for  Vx  which 
approaches  zero  as  x -*■  ± o°,  so  the  plasma  is  stable  in  ideal  MHD. 

Another  way  to  understand  this  ideal  MHD  stability  is  as  follows. 
According  to  Maxwells  equation  an  electric  field 

- 7TJ-  V,  (TOT  7) 

is  induced  by  the  fluid  motion.  However,  we  also  have 

± V*  Boj  (jUL  2 ) 

so  that  any  non  zero  E induces  an  infinitely  large  V near  the  singular 

z x 

point  x = = 0.  The  only  possible  resolution  in  ideal  MHD  is  of  course 

Y = 0,  or  stability  for  the  fluid.  If  the  resistivity  is  non  zero,  then 
near  the  singular  point  the  electric  field  can  be  balanced  not  only  by 
fluid  flow,  but  by  Ohmic  dissipation. 

Let  us  now  imagine  that  the  resistivity  is  very  small,  but  non- 
vanishing, and  that  its  presence  gives  rise  to  instabilities  whose  growth 
rate  vanishes  in  some  way  as  n -*■  0.  Then,  away  from  the  singular  point  we 
expect  the  plasma  to  be  described  by  the  Ideal  MHD  equations  with  zero  y, 
or  setting  the  right  hand  side  of  Eq.  (VII  4)  equal  to  zero. 
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Let  us  briefly  discuss  the  significance  of  the  neglect  of  y.  Neglecting 
the  growth  rate  in  the  momentum  equation  simply  means 

- (J  x B ) + - (J  x B)  - Vp  = 0, 
c - -o  CO 

or  equivantly 

V X | J X Bc  + Jo  X B | =0. 

That  is  away  from  the  singular  region,  the  perturbed  plasma  is  in  a state 

of  ideal  MHD  pressure  balance  equilibrium.  Imagine  that  the  tearing  mode 

we  consider  grows  with  growth  rate  much  less  than  kV  where  V is  a characteristic 

A A 

Alfven  speed  in  the  outer  region.  Then  if  the  plasma  configuration  is  not 
in  pressure  balance  equilibrium,  magnetosonic  waves  and/or  shocks  with 
frequency  oj  ^ kVA  will  be  generated.  Once  these  magnetosonic  disturbances 
damp  out,  pressure  balance  will  be  attained.  Our  fundamental  assumption 
then  is  that  the  time  scale  for  maintaining  this  pressure  balance  is  much 
less  than  the  characteristic  growth  time  of  the  tearing  mode.  Therefore, 
in  the  outer  region,  pressure  balance  applies  for  all  time  and  this  outer 
region  is  described  by  Eq.  VII  9. 
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There  are  two  difficulties  in  describing  the  entire  system  with  Eq.  (VII  9) 
above.  First  of  all  there  is  no  eigenvalue  y,  so  that  a solution  which 
satisfies  the  proper  boundary  condition  at  x = - <*>  in  general  will  not 
satisfy  the  boundary  condition  at  x = + Stated  another  way,  if  the 
solution  satisfies  the  boundary  conditions  at  both  x = - «>  and  x = + co^ 
in  general  there  will  be  a discontinuity  in  slope  at  the  position  where 
the  two  solutions  meet  (assuming  the  solutions  are  normalized  so  that 
there  is  no  discontinuity  in  itself).  Secondly,  the  equation  is 

singular  at  x = 0 where  B = 0.  As  shown  in  Appendix  , the  two  linearly 

independent  solutions,  denoted  ip  and  <p  behave  for  small  x as 

0$  = c+ 

The  idea  then  is  that  Eq.  (VII  1 ) applies  for  all  x except  a very  narrow 
range  around  x = 0,  (the  outer  region).  The  solution  of  Eq.  (VII  % ) in 
the  outer  region  has  a discontinuity  in  slope  at  x = 0.  In  a very  narrow 
range  around  x=  0 (the  inner  region),  inertia  and  resistivity  both  come 
into  play  and  smoothly  connect  the  two  outer  solutions. 

Before  actually  calculating  the  growth  rate  we  consider  more  carefully 
just  what  drives  instability  and  calculate  qualitative  expressions  for  the 
growth  rate.  Let  us  calculate  the  power  liberated  by  the  fluid  in  its 
motion.  Imagine  the  fluid  is  acted  upon  by  a force  density  F and  that  the 
fluid  velocity  is  V.  Then  the  power  P going  into  the  fluid  is 

v*-f  *c.c.  (jnr.  a ■) 

However  for  incompressible  motion, 

V = curl  R 


I'M  /p) 
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so  that  the  power  going  into  the  fluid  is 

p-  lOv  f (ySL  ii) 

where  we  have  assumed  the  solution  is  well  behaved  at  x = ± <*>  and  that  the 
velocities  at  x = 0+  are  somehow  smoothly  connected  to  those  at  x = 0 . 

For  two  dimensional  incompressible  motion 

^ nr-’ 


so  that 


P-  J'x 


-C  Jt 


-t  c.c 


W H 


In  ideal  MHD,  V is  related  to  B by  Eq.  (VII  6),  with  n = 0,  and  the 
x x 

right  hand  side  of  Eq.  (VII  4)  is  an  expression  for  the  z component  of 
the  curl  of  the  force  density.  Putting  these  together  we  find 

P-  - £[*  jr  & [-  + | J 

where  we  have  assumed  that  B is  real.  The  quantity  in  the  curly  brackets 

x 

is  just  Eq.  (VII  10)  and  is  zero  in  the  outer  region.  However,  Eq.  (VII  10) 
is  not  satisfied  everywhere  in  the  plasma,  but  only  in  the  outer  region. 

Thus  power  can  be  dissipated  by  the  fluid  in  the  inner  region. 

The  outer  region  solutions  for  large x generally  have  a discontinuous 
derivative  at  x = 0,  but  as  is  apparent  from  Eqs.  (VII  11  a and  b) , B 

x 

itself  is  well  behaved  and  continuous.  This  means  that  in  the  inner  region, 

g2g 

is  very  large  so  that  it  is  the  dominant  term  in  Eq.  (VII  5).  Hence 

"5x2 


the  power  per  unit  area  dissipated  in  the  inner  region  is  given  by 

p-  (zmo>  4 + r- 


(yn  a) 
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where 


Thus  if  A is  positive,  energy  is  dissipated  by  the  fluid  in  the  inner 

region,  and  we  expect  that  the  plasma  is  unstable.  Clearly  the  potential 

energy  driving  this  dissipation  comes  from  the  outer  region.  This  is 

obvious  because  B is  not  involved  in  the  motion  and  B is  zero  at  the 
oz  oy 

singular  point,  so  there  is  no  free  energy  in  the  inner  region.  In 
Chapters  VIII  and  XII  we  will  discuss  other  aspects  of  the  free  energy 
which  drives  the  tearing  mode. 

The  basic  picture  of  a tearing  mode  is  then  that  energy  is  released 
in  the  outer  region  and  dissipated  in  the  inner  region.  The  condition  for 
release  of  energy  in  the  outer  region,  and  therefore  the  condition  for 
instability  is 


A > 0. 


(VII  18) 


As  is  shown  in  Appendix  , for  the  sheet  current  at  x = ± a,  and  the 

prof ile 

B = B tanh  ka,  (VII  19) 

oy  oy  ' 


the  a's  are  given  respectively  by 

-Jka. 


4 - 


J?a  - (l- -fra.)  tannh  JpA 


(yjr 


IT 


Hence  for  sufficiently  small  k,  each  current  distribution  is  potentially 
unstable  to  tearing  modes. 
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One  other  interesting  aspect  of  the  instability  criterion  A > 0 
is  the  following.  If  A = 0,  the  MHD  solutions  for  y=0  connect  smoothly 
to  each  other  through  the  singularity.  Thus  if  A = 0 there  is  another 
MHD  equilibrium  (B  = Bq  + B)  which  can  be  arrived  at  from  the  original 
equilibrium  with  no  expenditure  of  energy.  Hence  the  condition  for  tearing 
mode  instability  is  that  there  be  a neighboring  equilibrium  of  the  magnetic 
structure  which  has  the  same  energy  as  the  initial  equilibrium.  This 
aspect  will  be  important  later  on  when  we  discuss  the  nonlinear  theory. 

Now  that  we  have  the  condition  for  instability,  we  continue  by 


estimating  the  size  of  the  inner  region  and  the  growth  rate.  To  start. 


we  will  eliminate  the  x term  in  Eq.  (VII  4)  by  using  Eq.  (VII  5). 
The  result  is 

V ' V/  lr‘  ***  j 


(yjr  sz) 


X 

where  we  have  neglected  k2B  <<  _ in  Eq.  (VII  5).  Within  the  inner 

x ox*- 

region,  B^  is  nearly  constant  as  we  have  discussed  previously,  and 

BQy  ~ Bq  - — . Thus  Eq.(VII  22)  is  an  inhomogeneous  differential  equation 
s 

for  Vx  in  terms  of  the  now  constant  B . Clearly  the  characteristic  length 


To  derive  the  growth  rate,  let  us  assume  that  the  energy  released 

in  the  outer  region  is  dissipated  by  Ohmic  heating  in  the  inner  region. 
3B 

Recall  that  x is  discontinuous  across  the  inner  region.  Since 

3x 
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1 3B 

By  = - , By  is  then  itself  discontinuous  across  this  region.  Thus 


there  is  a strong  current  in  the  z direction  which  flows  in  the  inner  region. 


A rough  estimate  of  the  magnitude  of  this  current  is 


C 

4-rf 


b!-  £ 


c 


l, 


, 7 


(yw  «,v) 


The  power  dissipated  per  unit  area  is  then  n Lc . Equating  the  Ohmic  power 


dissipated  in  the  inner  region  with  the  power  released  by  the  fluid  in  the 

y , 


outer  region,  we  find  z 


ihrJk' 


(jHL 


or 


and 


yz 


(jl 


(jhl~  cps-c.  j 


Notice  that  once  again,  as  for  the  case  of  the  resistive  g mode,  the  growth 
rate  is  proportional  to  a fractional  power  of  the  resistivity.  This  means 
that  even  if  7 x q V x Bq  ^ 0,  so  that  the  plasma  is  not  initially  in 
equilibrium,  the  time  for  the  relaxation  of  the  plasma  due  to  tearing  modes, 
t a,  n~3/5  is  quicker  than  the  time  scale  for  relaxation  of  the  assumed 
equilibrium  due  to  resistive  diffusion  t n~*. 

We  now  conclude  this  chapter  by  actually  calculating  the  growth  rate 
and  also  Ohmic  and  kinetic  power  dissipated  in  the  inner  region.  Here  the 
velocity  is  the  solution  of  Eq.  (VII  22).  We  now  make  a slight  simplification 
and  assume  that  the  current  sheet  is  symmetric  so  that  B " also  vanishes  at 
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the  singular  surface.  In  this  case,  in  the  limit  of  n 4 0 the  dominant 


g 

term  on  the  right  handside  of  Eq.  (VII  22)  is  the  - term  since  it 

nc 

diverges  as  n-1  as  n -*■  0,  whereas  the  other  two  terms  diverge  as  n-3/5 
according  to  Eq.  (VII  25b). 

The  idea  now  is  to  solve  Eq.  (VII  22)  for  V^  in  the  inner  region. 
The  growth  rate  is  then  calculated  by  insisting  that  the  inner  region 
solution  matches  with  the  outer  region  solutions  for  large  x.  If  the 
inner  region  is  narrow  on  the  length  scale  which  B varies  on,  then  B 

X X 

on  the  right  hand  side  of  Eq.(VII  ) can  be  regarded  as  constant.  This 
is  equivalent  to  the  constant  approximation  first  introduced  by  Furth, 
K ileen  and  Rosenbluth  (Physics  of  Fluids,  6^,  459  (1963)). 


Making  these  approximations,  Eq.  (VII  22)  reduces  to 

b>.v,-i^v,=  u (w-h 

' ^ I * fc re*  t-s 

assuming  Pq  does  not  vary  substantially  in  the  inner  region.  The  solution 

for  V^  from  Eq.  (VII  26)  is  that  linear  combination  of  homogeneous  solutions 

and  the  particular  solutions  which  satisfy  the  appropriate  boundary 

conditions,  that  is  V -*•  0 as  x -*■  ± ».  We  first  look  at  the  solution  to 

x 

the  homogeneous  equation 


Al\/  = / a t ) V 

Put  ( i 


Since  the  quantity  in  parentheses  is  positive,  the  solution  for  V^_  which  is 
positive  at  x + - <•  is  a monotonically  increasing  function  of  x and  there- 
fore cannot  possibly  satisfy  the  boundary  conditions.  Therefore  the 
appropriate  solution  to  Eq.  (VII  26)  is  the  particular  solution  alone. 

In  general  this  can  easily  be  calculated  numerically.  Here  we  give  an 
approximate  analytic  solution.  In  the  limit  of  q -*■  0,  the  particular 
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solution  for  large  x is 

t / Li  5* 


t x 


V*  * 


A K* 


j-  a 


a 


('SEji') 


tf 


which  is  the  MHD  solution.  We  now  proceed  to  calculate  the  particular 
solution  for  all  x.  Assuming  that  the  inner  region  width  is  much  less 
than  1/k,  Eq.  (VII  ) can  be  simplified  to 


L, 


r dx(^J  U * 


A So  Lc 


At  large  x,  the  particular  solution  is  given  by  Eq.  (VII  27)  while  for 
small  x it  is  the  solution  to 


J / 


c y 3X  (*=*)  L- 

~A  30  iT 


Let  us  assume  a solution  for  Vx  of  the  form 


Ax 


/+  Cxl 

For  large  x,  this  matches  with  Eq.  (VII  27)  if 

£_  _ _ 8,  Jj) 

c J?, 


(jut  s.) 


(jnc  31 ) 


Expanding  Eq.  (VII  30)  for  small  x and  taking  the  second  derivative,  we 
find  that  Eq.  (VII  30)  matches  the  solution  for  small  x if 

Y 


Ot  3±) 


£ A 5 Cc  * 

Solving  for  A and  C,  we  find  an  approximate  particular  solution  to 


Eq.  (VII  28)  which  matches  the  approximate  solution  for  both  large  x and 
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small  x is 


yjj_  x 15x(x*°) 

J 6 A (yrr  33 ) 

i r Lt 

Thus  we  now  have  an  approximate  solution  for  Vx  in  the  inner  region. 

The  final  task  is  to  use  Eq.  (VII  5)  to  compute  the  eigenvalue  y.  The 
trick  is  to  insure  that  the  A generated  by  the  inner  region  solution  is 
the  same  A generated  by  the  outer  region  solution,  for  instance  Eqs.  (VII 
20  or  21).  To  calculate  the  delta  generated  by  the  inner  solution,  we 
utilize  Eq.  (VII  5 ) (with  k2Bx  neglected  on  the  right  hand  side).  Inserting 
Vx  from  Eq.  (VII  33),  this  is 

dll*s  *_r 1 

J*  fc  L '+  & J 

Integrating  Eq.(VII  34)  from  x=-°°tox  = + <*>  (that  is  across  the  inner 
region),  we  find 

. &)\  J&- (£)*  *£* 

* ' -co 

Inserting  for  L from  Eq.  (VII  23),  we  see  that  the  growth  rate  is  given 

to  within  numerical  factors  by  the  approximate  expression,  Eq.  (VII  25b). 

We  now  calculate  the  Ohmic  power  dissipated 

in  the  singular  layer,  and  the  kinetic  power  dissipated  there  also.  The 

Ohmic  power  dissipated  per  unit  area  is  qJ  2 integrated  across  the  singular 

z 

layer.  Now  since 

£z  - (tt3i) 
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but  is  also  equal  to  -yB  /ike;  we  find  that  the  Ohmic  power  dissipated  is 


The  power  dissipated  into  kinetic  energy  is  somewhat  more  difficult 
to  define  because  there  is  kinetic  power  in  both  the  inner  and  outer  region. 
However,  as  is  apparent  from  Fig.flll  1)  and  also  from  Eq. (VII  3),  the  y 


velocity  is  much  larger  than  the  x velocity  in  the  inner  region,  but  much 
smaller  than  it  in  the  outer  region.  Therefore,  we  will  define  |pV^  as  the 
kinetic  energy  dissipated  in  the  layer.  Using  the  fact  that  V*V  = 0,  the 


kinetic  power  dissipated  in  the  layer  is  then  given  by 


where  we  have  made  use  of  Eq.  (VII  33)  for  V . The  integrals  in  Eqs.  (VII 
38  and  39)  are  forms  (or  slight  variations  of  forms)  which  can  be  found  in 
tables  of  definite  integrals  for  instance 

Using  the  expression  for  y and  L , we  find 

Tk  j_ 

/( 
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I Therefore,  in  a tearing  mode,  most  of  the  energy  liberated  by  the 

fluid  in  the  outer  region  is  dissipated  into  Ohmic  heating  of  the  inner 
region,  with  very  little  going  into  fluid  motion  in  the  inner  region, 
u To  summarize,  the  tearing  mode  just  becomes  unstable  when  there  is 

a second  magnetic  equilibrium  which  is  accessible  from  the  first  with  no 
expenditure  of  energy. 


The  instability  develops  by  the  fluid  releasing  magnetic  energy  in  the 
outer  region,  and  dissipating  it  principally  via  Ohmic  heating  in  a narrow 
inner  region  near  the  singular  surface.  This  dissipation  near  x = 0 manifests 
itself  by  a change  in  the  magnetic  field  topology,  that  is,  island  formation. 

If  By<<BQ,  the  equations  for  the  field  lines  in  the  xy  plane  is 




B cos 


3, 


'*  S-s 


where  Bx  is  taken  as  Bcosky.  The  equation  for  the  field  line  is  then 


((VII  40) 


2L  B 

s 

The  field  line  has  a separatrix  if  K = and  this  separatrix 

o 

has  maximum  width  (the  island  width) 


(VII  41) 


The  field  lines  in  the  xy  plane  are  shown  in  Fig.  (VII  4). 
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Fig.  VII 


89 


90 


VIII.  Internal  MHD  Instabilities  In  Cylindrical  Plasmas 


In  this  section  we  discuss  MHD  instabilities  in  a cylindrical  plasma 
with  no  free  surface.  As  we  will  see  shortly,  the  modes  which  we  have 
discussed  in  the  previous  three  chapters  have  simple  analogs  in  cylindrical 
geometry.  However  there  are  also  additional  modes  in  cylindrical  geometry 
which  have  no  analog  in  slab  geometry.  Because  this  chapter  is  quite  long, 
we  subdivide  it  into  six  subchapters,  VIII.  A - VIII.  F,  which  deal  respectively 
with  two  dimensional  MHD  modes,  ideal  MHD  modes  with  m » 2,  ideal  MHD  modes 
with  m = 1,  resistive  MHD  modes  with  m > 2,  resistive  MHD  modes  with  m = 1, 
and  double  tearing  modes. 

VIII.  A - MHD  Modes  with  Two  Dimensional  Structure 

Recall  from  the  last  three  chapters,  that  in  slab  geometry,  the 
appropriate  modes  always  had  two  dimensional  structure  in  the  pla*.e 
perpendicular  to  B . If  Bz0>>Bg0  and  a<<R,  (tokamak  ordering)  one  might 
expect  this  to  be  valid  in  cylindrical  geometry  also.  This  is  analogous  to 
the  analysis  in  the  second  part  of  Chapter  IV  tfiere  we  examined  two  dimensional 
unstable  modes  in  a cylindrical  plasma  with  a free  surface.  Here,  as  we  will 
see,  two  dimensional  modes  in  a cylindrical  plasma  with  no  free  surface  are 
stable.  However  the  analysis  here  is  still  useful  to  set  the  stage  for  a 
study  of  instabilities  resulting  from  three  dimensional  or  resistive  effects. 

In  cylindrical  geometry,  all  perturbed  quantities  vary  as 
exp  i(m9  + kz).  To  make  the  transition  to  toroidal  geometry  we  simply 
quantize  k,  that  is  k + n/R,  but  neglect  all  other  complications  of  the 
toroidal  geometry.  The  perturbation  is  perpendicular  to  the  ambient  field 


wherever  mB^r)  + nB  =0,  so  we  expect  unstable  perturbations  to  be 

X ZX 
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in  some  way  localized  near  this  point  (we  assume  here  that  B is  independent 

B Z° 

of  r and  is  >>B„  ).  If  this  is  so,  k « ^ ^ <<  ^ so  that  the 

variation  in  the  9 direction  is  much  more  rapid  than  the  variation  in  the 
z direction.  Therefore,  plasma  motion  will  be  approximately  2 dimensional 
in  r and  9,  quite  analogous  to  the  motion  in  slab  geometry.  Then  = 0, 

k is  taken  as  a small  quantity  unless  it  multiplies  B^  and  V^,  V^,  and 

B,  are  related  by  the  incompressibility  condition 

y 


/ 9.  r v/  + c y*t  . / , _/  9 r p j-  LlS:  V.  - C> 

r 2>r  > T ' r br  r r ® 


(VIII  1) 


The  next  step  then  (analogous  to  the  previous  chapters)  is  to  take  the  z 
component  of  the  curl  of  the  perturbed  momentum  equation,  or 

j(8-  vJr.  + (B. • ?j  X - (?•  ?)  ?.  - (T.  • 5)  ? j 


_L 

<£ 


(VIII  2) 


First  look  at  the  third  and  fourth  terms  on  the  right  hand  side  of  Eq.  (VIII. 2) 

Since  B is  assumed  to  be  constant,  (J*V)B  = 0.  The  fourth  term  is  - ikJ  B 
oz  o o z 

which  is  small  since  both  k and  B are  assumed  small.  The  first  term  on 

1 dJoz  Z im 

the  right  is  — B while  the  second  is  ( — B„  + ikB  )J  . However 

0 c r dr  r 0o  zo  z 


J - f- 
z 4ir 


VxB^  . Making  use  of  Eq.  (VIII.  1),  Eq.  (VIII.  2)  reduces  to 


I i *7  L£  cL  i2. 


rdr'^Zr 


(VIII  3) 


where 


F= 


yrt 


'&o 


(VIII  4) 


Here,  unlike  the  second  part  of  Chapter  IV,  F depends  on  r in  general. 
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To  complete  the  description  of  linear  modes  in  a cylindrical  plasma  in 
the  tokamak  ordering,  an  additional  relation  between  V^_  and  is  required. 
This  comes  from  Maxwell's  equation 


Y 8 - - v * y * B 


V*  ^ * 3 


(VIII  5) 


To  start  we  consider  only  ideal  modes,  so  q = 0 in  Eq.  (VIII.  5) 
above.  In  this  case 

y 7?  » i F vr 

0 r (VIII  6) 

Then  inserting  for  on  the  right  hand  side  of  Eq.  (VIII.  3),  we  find 


after  a bit  of  straight  forward  manipulation 

J-  jL 
r 


r’/~  r r^Vr 

1 (VIII  7) 


Equation  (VIII.  7)  above  must  be  solved  subject  to  the  boundary  condition 

that  V = 0 at  r = a and  V is  well  behaved  for  r = 0. 
r r 


Near  r = 0,  Eq.(VIII.  7)  reduces  to 

= - r(l-  ^)vr 

?r 


(VIII  8) 


so  that  the  solutions  have  the  form  V^  = r . Solving  for  y*.  (See  Appendix 


), 


we  f ind 


- -)  ± 


(VIII  9) 


For  m > 2,  one  root  is  well  behaved  at  r = 0 and  one  is  not.  The  solution 
which  is  well  behaved  at  zero  goes  as  r1  ^ so  that  V = 0 there.  Hence  t • 
center  of  the  plasma  does  not  move. 

On  the  other  hand,  if  m = 1,  the  well  behaved  solution  f 
constant  at  r = 0.  Let  us  now  see  what  this  implies  for  the  n. 
plasmacenter . If  4r  is  a unit  vector  in  the  r direction  V^  = 
r = 0,  V^_  is  a constant  times  cos(9  + kz)  and  j_r  = L cos  * 

Therefore  at  z = 0,  the  velocity  V is  in  the  positive  • 
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for  instance  z = the  velocity  V is  in  the  negative  y direction.  Clearly 
then,  an  m = 1 mode  corresponds  to  a helical  displacement  of  the  center  of 
the  plasma. 

Since  modes  with  m » 2 do  not  displace  the  center  of  the  plasma, 
the  cylindrical  geometry  is  not  crucial  and  modes  with  m » 2 are  not  very 
different  from  analogous  modes  in  slab  geometry.  On  the  other  hand  modes 
with  m = 1 represent  a rigid  helical  displacement  of  the  center  of  the  plasma. 
Since  there  is  no  analogous  motion  in  a slab,  we  find  that  m = 1 modes  in  a 
cylinder  can  be  very  different  from  anything  occuring  in  a slab.  We  now 
examine  whether  Eq.  (VIII.  7)  predicts  instability. 

If  m >2,  the  coefficient  of  V^.  on  the  right  hand  side  of 
Eq.  (VIII.  7)  is  positive  for  p'<0,  so  that  if  V^  is  a monotonically 
increasing  function  of  r as  long  as  it  is  positive  in  the  viscinity  of 
r = 0.  Thus  there  is  no  solution  to  Eq.  (VIII.  7)  which  satisfies  the 
boundary  condition  at  the  wall,  V^_(r  = a)  = 0;  hence  there  are  no  unstable 
modes  with  m > 2. 

If  m = 1,  the  solution  to  Eq.  (VIII.  7)  which  is  well  behaved  at 
Che  origin  is  V^_  is  constant.  Since  V^_  must  vanish  at  the  wall,  it  vanishes 
everywhere,  so  there  are  no  unstable  modes  with  m = 1 either. 

Hence  the  conclusion  is  that  in  a cylinder,  there  are  no  unstable 
modes  in  ideal  MHD  which  have  two  dimensional  structure  and  no  free  surface. 

The  reason  is  that  the  coefficient  of  V^  on  the  right  hand  side  of  Eq.  (VIII.  7) 
is  always  non  negative,  so  that,  as  discussed  in  Chapter  V,  there  can  be  no 
instability.  However  while  this  coefficient  is  not  negative,  it  can  be 
nearly  zero  wherever  F = 0,  that  is  wherever  the  perturbation  is  perpendicular 
to  the  ambient  field.  Thus  small  corrections  to  either  the  two  dimensional 
structure  of  the  mode,  or  to  ideal  MHD  (that  is  resistive  effects)  might  give 
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rise  to  instabilities  which  are  somehow  centered  near  the  positions  where 
F = 0.  In  the  next  five  subsections,  we  will  see  that  this  is  indeed  the 
case. 


VIII.  B - Ideal  MHD  Modes  with  m » 2 


In  this  subsection  we  consider  ideal  MHD  modes  with  m > 2.  We 
deal  with  the  full  three  dimensional  mode  structure  and  make  no  approximation 
concerning  either  Bg/B^  of  kr/m.  The  starting  point  is  the  linearized  fluid 


equation  and  Maxwell's  equation: 


' ’ 9*  J V-rr-  </r 

r 8 -■  - vy 


fyrrr  i° 

(ynr.  * i 


for  the  two  vector  quantities  V and  B.  Consistent  with  our  usual  notation, 
a quantity  without  a subscript  is  a perturbed  quantity  and  a quantity  with 
a zero  subscript  is  an  equilibrium  quantity.  The  perturbed  pressure  is 
eliminated  by  the  incompressibility  condition 


The  idea  now  is  to  reduce  this  set  of  equations  to  a single  equation  for  V . 
The  manipulations  involved  in  this  are  tedious  and  are  set  out  in  Appendix  _ 
The  final  result  is 


where 


^ /»  etVr  __  by r -o 
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(jmc  i?c ) 


and  where  F is  given  by  Eq.  (VIII.  4). 

Notice  that  as  long  as  y2>0,  Eq.  (VIII.  12  a)  is  nonsingular  and 

one  could  simply  solve  it  numerically  for  y.  This  is  the  approach  taken 

by  Freidberg  (Physics  of  Fluids  13,  1812  (1970)).  However  there  are  many 

insights  that  can  easily  be  obtained  analytically.  Specifically,  notice 

that  Eq.  (VIII  12  a)  is  very  similar  in  structure  to  Eq.  (V.  17),  derived 

in  our  discussion  of  gravitational  modes  in  a sheared  field.  That  is, 

regions  of  negative  h are  unstable  and  regions  of  positive  h are  stabilizing. 

As  long  as  < 0,  h will  always  be  regative  (for  small  y)  at  least  at  the  singular 

m2B§ 

surface  where  F ■ 0 (of  course  k2B2 y—  = 0 where  F - 0) . 

z r4. 

As  in  the  case  of  g modes  in  a slab,  for  fixed  y(=yQ)  there  will 
always  be  an  unstable  mode  with  y > yQ  as  long  as  the  eigenfunction  V^(r) 
has  a node  between  r * 0 and  r * a.  The  first  possibility  then  is  that 
Vf(r)  has  a node  near  F ■ 0 in  the  limit  of  y -+■  0.  A calculation  identical 
to  that  leading  up  to  Eq.  (V.  30)  shows  that  such  a localized  mode  will  be 
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rB 

where  we  have  used  the  fact  that  q = — — and  k = for  a toroidal  plasma. 

RB  K 

9o 

Also,  at  the  singular  surface  q = - This  is  the  Suydam  condition  for 

instability  (Suydam,  Geneva  UN  Atomic  Energy  Conference,  1958  p.  1057). 

A large  q'  can  stabilize  this  mode,  as  shown  in  Eq.  (VIII.  13). 

To  see  that  this  is  a shear  stabilization,  let  us  calculate  the  shear  length 


in  cylindrical  geometry.  Recall  that  in  slab  geometry,  k,,  = k^  7^.  In 

zo  11 


cylindrical  geometry,  then  k = 


(Bo*7)  \ B0o  + % BZ 
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Near  the 


position  r where  F = o,  kii  = [ — p— 5 — - —7—77—  — -j^-1  (r  - r ) 
r s ’ | r [B2  + bz  q dr  s 
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so  that  if  ^ is  taken  to  correspond  to  k^,  then 
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Now  let  us  relate  these  pressure  driven  modes  to  the  gravitational 
modes  discussed  in  Chapter  V.  Setting  = F = 0 in  Eq.  (VIII.  12  a)  and 


taking  the  limit  m2>>k2r2,  it  reduces  to 

2 )*V 

% Por  + \ ^rrr  / _ 
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Clearly,  a negative  pressure  gradient  causes  an  instability  whose  growth  rate 
is  roughly 


?:  Y & 
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However  since  F » 0,  k£/  ■ • The  radius  of  curvature  of  the  field 

m ' zo 

line  is  given  by  2. 

j_  _ . (3.-v)3.  _ j_  3».  i 
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and  is  directed  outward.  Thus  the  unstable  mode  is  like  a gravitational 

mode  with  g yp-  given  roughly  by  — ^ where  R is  the  radius  of 

M o 

curvature  of  the  field  line.  This  is  very  similar  to  that  derived  in 
Eq.  (V.  5). 

To  summarize,  we  have  shown  that  modes  very  similar  to 
gravitational  modes  can  be  driven  unstable  by  a negative  pressure  gradient. 
These  can  be  stabilized  by  shear  as  shown  in  Eq.  (VIII.  13).  However 
the  plasma  is  not  necessarially  stable  if  Eq.  (VIII  13)  is  violated  at 
every  point.  Equation  (VIII  13  ) is  only  the  condition  for  unstable  modes 
localized  about  a particular  singular  surface.  There  also  may  be  gross 
modes  which  are  not  localized.  This  possibility  has  been  investigated  by 
Newcomb  (Annals  of  Physics  10,  232  (I960)).  He  examined  when  the  solution 
of  Eq.  (VIII  12  a)  can  have  a node  in  the  limit  of  small  y.  One  can  of 
course  also  simply  solve  Eq.  (VIII  12  a)  numerically,  the  approach  taken 
by  Freidberg. 

VIII.  C - Ideal  MHD  Modes  with  m - 1 

For  a cylindrical  plasma  in  Heal  MHD,  one  can  always  determine 

mode  structure  and  growth  rates  by  solving  Eq.  (VIII.  12  a)  subject  to 

proper  boundary  conditions  at  r = 0 and  r * a.  If  BQo  ^ Bzq  and  kr  ^ m, 

as  in  for  instance  a reversed  field  pinch,  there  is  no  particular  distinction 

between  modes  with  m * 1 and  modes  with  all  other  m.  However  in  tokamak 

ordering,  BQ  <<B  , kr<<m,  the  behavior  of  modes  with  m - 1 is  strikingly 

oo  zo 

different  from  the  behavior  of  modes  with  all  other  m. 

Recall  from  Chapter  V that  only  those  regions  of  the  plasma 
having  h < 0 (h  is  defined  in  Eq.  (VIII.  12  c)  can  drive  instability.  In 
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tokamak  ordering  for  m ^ 2,  the  second  term  in  h is  so  large  that  h can 
only  be  negative  in  a small  region  around  the  singular  surface.  However 
for  m = 1,  this  term  nearly  vanishes.  If  k is  equal  to  then  in 

general,  for  m = 1 and  q(r)  a monotonically  increasing  function  of  r,  h 

will  be  negative  from  r =*  0 to  the  radial  position  at  which  q */.  Hence- 

7t 

forth  we  specialize  to  the  case  of  Tl  = 1»  so  that  the  unstable  region  is 
within  the  q = 1 surface. 

It  was  shown  in  Section  VIII.  A that  motions  with  m * 1 are 
unique  in  another  way,  namely  these  are  the  only  motions  which  displace 
the  plasma  center.  Therefore  these  modes  are  expected  to  displace  the 
entire  central  region  of  the  plasma.  We  continue  by  calculating  the  growth 
rates  and  eigenfunctions  for  these  modes  and  conclude  by  discussing 
physically  what  is  happening. 

In  the  limit  of  low  growth  rate,  we  neglect  the  y2  terms  in  h. 
Thus  the  equation  describing  the  motion  is  still  given  by  Eq.  (VIII  12a) 
with 


Notice  that  h is  multiplied  by  (✓£) 2 which  is  a small  term.  Therefore,  Vf 
is  given  approximately  by  the  solution  to  the  equation 

± r‘  U A "?)  J7 - ° CVZ  18  ) 

I 

As  we  will  see  shortly,  V^  is  nearly  constant  for  r within  the  singular 

dVr 

surface,  and  is  nearly  zero  for  r outside.  Thus  is  nearly  zero 
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everywhere  except  in  the  viscinity  of  the  singular  surface  r«rg. 


Anticipating  this  result,  and  setting  q-1 
Eq.  (vill  it) 

oUrr 


q ' (r-r  ) , we  find  from 
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Con  s 


(VIII  19) 


(VIII  20) 


and  the  two  different  values  are  smoothly  connected  in  a narrow  region  of 


width 


4nV2rs 

Bf  (r  )q'2 
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about  r = rg.  The  proper  boundary  condition  of 


course  is  (r«a)  = 0,  which  dictates  the  relation  between  Kj  and  K2 

K,  --  - % *. 


(VIII  21 1 


Also,  we  find 


K-,- 


- W Vr(r*o) 


(VIII  22) 


The  eigenfunction  V^(r)  is  shown  in  Fig.  VIII.  1. 

Now  we  must  solve  for  y.  Of  course  y must  result  from  the  h term 

which  drives  the  instability.  Assuming  that  Vf  is  constant  between  r = 0 

and  r » r , we  find  that  another  relation  for  dVr  is 
s — 3— 

dr 


oivr 
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Of  course  the  — from  Eq.  VIII.  23)  must  be  consistent  with  that  derived 
in  Eq.  (VIII  19  or  20)  in  the  limit  as  r -*■  rg.  This  allows  us  to  solve  for 
Y,  the  result  being  ^ 


(VIII  24) 


r 

Thus,  as  long  as  f h(r)  < 0,  the  mode  is  unstable, 
o 

We  now  discuss  briefly  the  physical  mechanism  of  this  instability. 

As  emphasized  many  times,  the  most  unstable  perturbations  are  always  those 
which  bend  the  field  lines  least.  In  slab  geometry,  with  boundary  conditions 
imposed  at  ± °°,  this  means  that  the  perturbation  is  somehow  localized  near 
the  singular  surface.  In  cylindrical  geometry  however  there  is  another 
possibility.  A rigid  displacement  of  a cylinder  certainly  does  not  bend 
any  field  lines.  The  problem  is  that  a rigid  displacement  cannot  satisfy 
the  boundary  condition  V^(r=a)  = 0.  The  question  then  is,  can  only  the  inner 
part  of  the  cylinder  be  rigidly  displaced  but  still  not  bend  the  field  lines? 

To  investigate  this  question,  we  first  examine  the  fluid  motion. 

The  motion  V = constant  for  r < rg  is  of  course  something  of  an  over- 
simplification. Clearly  such  motion  cannot  be  incompressible.  To  describe 
more  precisely  this  incompressible  motion,  imagine  that  in  a narrow  layer 
of  width  5 about  r = rg,  V^.  decreases  linearly  to  zero.  In  this  case 


o,  r<rs  J v, 


-t 


S <r<rs+S  j 
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rs+S-r 


- Vr 


(VIII  25) 
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This  motion  is  reminiscent  of  the  tearing  motion  described  in  the  last 

chapter.  Assuming  that  the  radial  velocity  goes  as  cos(m0  + — ) , the 

R 

velocity  streamlines  are  shown  in  Fig.  (VIII  2)  at  four  different  axial 

ttR  3^ 

positions,  Z = 0,  ^ , ttR  and  • The  streamlines  are  solid  lines, 

while  the  two  circles  r = r and  r = r +6  are  dashed.  Also  shown  on 

s s 

each  plot  are  three  dots  labled  1,  2 and  3.  Dot  number  one  corresponds 

to  an  initial  position  just  inside  the  singular  layer,  where  q = 1.  Since  \ 

the  pitch  of  the  field  line  is  the  same  as  the  pitch  of  the  perturbation, 
the  same  field  line  passes  through  all  dots  labeled  1.  As  the  flow  proceeds, 
for  a time  and  t.  y the  fluid  element  passes  to  positions  2 and  3,  while 
the  field  line,  of  course,  continues  to  interlace  these  points.  In  each  of 
the  four  curves  shown,  the  displacement  between  positions  1 and  2 and 
between  2 and  3 is  the  same  except  for  a 90°  rotation.  Thus  the  field  line 
is  not  bent  or  stretched  as  the  fluid  races  around,  from  front  to  back, 
near  the  singular  surface.  For  r < rg  the  fluid  motion  is  a rigid  displacement 
in  the  (r9)  plane  which  does  not  bend  field  lines.  For  rg<r<rg+<5,  the 
pitch  of  the  field  line  is  the  same  as  the  pitch  of  the  perturbation, 
allowing  the  field  line  to  be  rigidly  displaced  around  the  singular  layer. 

Therefore  the  type  of  motion  shown  in  Fig.  (VIII  2)  does  not  initially  bend 
field  lines  and  is  likely  to  be  unstable. 

Now  examine  the  motion  for  a field  line  which  at£  = 0 is  not  near 
the  singular  surface.  In  the  standard  tokamak  configuration,  q < 1 for 
r < rg,  so  this  field  line  winds  a tighter  helix  at  £ = 0.  The  cylinder 
which  the  field  line  is  initially  on  is  dotted  in  Fig.  VIII  2 and  at 
the  various  Z’s  the  field  line  passes  through  the  points  A.  However,  the 
field  line  does  not  enter  the  singular  region  at  the  same  time  any  more. 
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since  at  some  axial  points  it  is  initially  closer,  along  a stream  line. 


to  the  singular  region.  Thus  some  points  of  the  field  line  whip  around 
the  singular  region  before  other  points  do.  Clearly,  this  involves  consid- 
erable bending  and  stretching  of  the  field  line  and  is  a stabilizing  effect. 
Therefore  as  the  motion  proceeds,  there  is  an  opposing  force,  which 
increases  with  the  displacement.  rather  like  compressing  a spring. 

The  distortions  of  a field  line  with q = 1 and  q < 1 are  shown  in  Figs. 

(VIII  3 and  4).  This  lets  us  interpret  somewhat  the  various  effects  which 
contribute  to  y in  Eq.  (VIII  24).  The  integral  h simply  represents  the 
average  force  which  displaces  the  plasma,  while  the  q'  in  the  denominator 
represents  the  restraining  force  from  the  bending  of  the  field  lines.  Of 
course  if  q " = 0,  the  analysis  leading  up  to  Eq.  (VIII  24)  is  invalid. 

For  q " = 0,  however,  the  analysis  leading  up  to  Eq.  (VIII  16)  is  still 
valid  and  we  find  that  for  zero  shear,  the  growth  rate  is  given  by 


(VIII  26) 


so  that  there  is  instability  if  Pq  and  pA  have  the  same  sign.  The  growth 
rate  in  Eq.  (VIII.  26)  is  much  greater  than  that  in  Eq.  (VIII.  24),  so  the 
presence  of  shear  slows  down  these  modes  hit  does  not  stabilize  them. 

VIII  D - Resistive  Instabilities  for  m > 2 

In  this  section  we  discuss  resistive  instabilities  in  a plasma 
with  m > 2.  The  two  types  of  instabilities  are  resistive  g modes  and 
tearing  modes.  In  section  VIII  b,  it  was  shown  that  in  a cylindrical 


3p 


o 1 


dP 


plasma,  there  are  simple  analogues  to  g modes  where  g ^ — — . Here 

o 

R is  the  radius  of  curvature  of  the  field  lines  and  the  sign  is  such  that 
a negative  pressure  gradient  gives  rise  to  instability.  In  ideal  MHD, 
these  pressure  driven  modes  could  be  stabilized  by  sufficient  shear,  as 
indicated  in  Eq.  (VIII  13).  However,  in  slab  geometry,  it  was  shown  in 


A 
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Chapter  VI  that  finite  resistivity  could  destabilize  these  modes  by 
allowing  the  plasma  to  leak  through  the  sheared  field.  Clearly,  in 


cylindrical  geometry,  we  expect  the  same  result  to  apply,  namely  that, 
analogous  to  Eq.  (VI  22),  finite  resistivity  allows  pressure  driven 
modes  to  go  unstable  with  a growth  rate  given  by 


H s 


l ^ \ r & 


(VIII  27) 


where  is  given  in  Eq.  (VIII  17)  and  Lg  in  Eq.  (VIII  14). 

The  tearing  mode  for  m > 2 in  tokamak  ordering  (B.  <<B  , kr<<l) 

bo  zo 

also  follows  in  a straight  forward  way  from  analogous  results  in  slab 
geometry  discussed  in  Section  VII.  That  is  if  the  inertial  terms  in 
Eq.  (VIII  3)  are  set  equal  to  zero,  the  equation  for  the  perturbed  field 
gives 


m2-  r-^r  5 - 


2TL  + 31  z 


3 , 


(VIII  28) 
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if  k = the  quantity  F = / (m-nq)  vanishes  at  the  singular  surface 

k r 

r = rg  where  q = %.  The  procedure  now  is  exactly  as  it  was  for  the  slab 

(Chapter  VII).  First  calculate  solutions  for  Eq.  (VIII  28)  for  r < rg 

and  r > rg  which  satisfy  the  proper  boundary  conditions  at  r = 0 and  r = a. 

If  these  solutions  are  normalized  to  have  the  same  value  of  B at  r = r , 

r s 

in  general  there  will  be  a discontinuity  in  slope  as  one  crosses  r . Then 
if  A (defined  in  Eq.  (VII  17))  is  positive  the  plasma  is  unstable  to  tearing 
modes.  Since  the  width  of  the  tearing  layer  is  very  small  compared  to  rg 
the  connection  between  the  two  MHD  solutions  can  proceed  exactly  as  in  slab 
geometry.  Thus  to  see  whether  a cylindrical  equilibrium  in  the  tokamak 
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ordering, is  stable  to  tearing  modes,  one  simply  solves  Eq.  VIII  28  numerically 
and  investigates  the  condition  for  A>0.  Calculations  such  as  these  have 
been  made  by  a variety  of  authors  (For  instance  F with,  Rutherford  and 
Selberg,  Phys.  Fluids  _16,  1054  (1973)). 

Stated  very  qualitatively,  the  conclusion  is  that  for  sufficiently  smooth 

current  profiles,  tearing  modes  are  £able  if  m >,  4.  For  m = 2 or  3,  the 

profile  is  typically  unstable.  However  by  tayloring  the  current  profile, 

it  is  possible  to  stabilize  these  modes  also. 

We  conclude  this  subsection  with  a simple  physical  picture 

illustrating  the  free  energy  which  drives  a tearing  mode  in  an  imcompressible 

iV*r 

plasma.  To  begin,  take  Eq.(VIII  3),  multiply  by  — and  integrate 
over  the  plasma  volume.  Integrating  the  second  derivative  terms  by  parts 


and  neglecting  the  end  point  contributions,  we  find  the  result: 


(VIII  29) 


In  calculating  Eq.  (VIII  29)  we  have  made  use  of  the  incompressibility 

condition  to  relate  V.  and  B.  to  V and  B . Notice  that  the  left  hand 

9 9 r r 

side  of  Eq.  (VIII  29)  is  simply  the  rate  of  change  of  kinetic  plus  magnetic 
energy.  The  term  on  the  right  hand  side  then  represents  a driving  term. 
This  term  can  only  drive  instability  in  those  regions  of  the  plasma  where 
dJo  1 dJ 

~ ^ and  F have  opposite  signs.  If  — — is  everywhere  positive,  then 
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one  can  easily  show  from  Eq.  (VIII  28),  that  for  boundary  conditions 

= 0 at  r = 0 and  r = a,  B^  is  a monotonically  increasing  function  of  r 
for  r < rg  and  a monotonically  decreasing  function  of  r for  r > rg. 

Thus  A<0  and  the  plasma  is  stable.  Hence  A">0  can  only  arise  if 
^ < 0,  somewhere  in  the  plasma. 

In  the  standard  tokamak  configuration  4^-  is  negative  and  q(r) 

r B9 

is  a monotonically  increasing  function  of  r.  Therefore  since  F = (m-nq), 

these  terms  have  opposite  signs  for  m > nq,  or  for  radii  less  than  the 
radius  of  the  singular  surface.  Thus  for  a tearing  mode  in  cylinder,  the 
region  inside  the  singular  surface  releases  free  energy;  the  region 
outside  soaks  some  of  it  up.  What  is  left  over  is  deposited  as  Ohmic 
hearing  near  the  singular  layer. 

It  is  instructive  to  examine  physically  the  nature  of  this  driving 
energy.  Since  the  plasma  motion  is  incompressible,  the  motion  of  each 
fluid  element  can  be  expressed  as  a displacement  plus  a rotation.  We 
will  consider  for  now  the  rotational  part  of  the  motion.  The  power  input 
into  rotation  is  the  torque  times  the  angular  velocity.  In  Fig.  (VIII  5) 
is  shown  a fluid  element  in  polar  co-ordinates  between  r and  r + dr  and 
9 and  9 + d9.  To  calculate  the  torque,  we  calculate  the  difference  between 
the  9 componentsof  the  force  on  the  two  edges  at  r and  r + dr.  Clearly  one 
element  of  this  torque  t = i x Force  is 


(VIII  30) 


The  problem  now  is  to  calculate  the  angular  velocity  of  the  fluid  element 
about  its  own  axis.  To  calculate  the  angular  velocity  about  the  axis,  first 
subtract  out  the  displacement  of  the  axis.  Then  an}  variation  of  radial 
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velocity  with  polar  angle  9 must  be  a rigid  rotation  about  the  axis.  Thus 


this  angular  velocity,  ft  is  simply 

C 2 o1  Vr 


oie 
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where  negative  — — 

ou 

relation  between  B 
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implies  fi  is  in  the  positive  Z. 
and  V in  ideal  MHD  is 


In  real  notation,  the 
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Therefore 
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so  that  F and  having  opposite  sign  means  that  P is  positive.  Now 

dr 

the  interpretation  is  clear  ;when  F and  dJoz/dr  have  opposite  signs,  the 
torque  on  a fluid  element  is  in  the  same  direction  as  the  rotation,  thereby 
tending  to  increase  the  rotation. 


VIII  E - Internal  m = 1 Tearing  Modes 

As  we  have  seen  in  Chapter  VIII  C for  tokamak  ordering,  modes  with 

m = 1 can  be  driven  unstable  in  ideal  MHD  by  the  three  dimensional  nature 

of  the  motion.  However,  as  noted  in  VIII  c the  basic  motion  is  two 

dimensional  and  it  consists  of  a rigid  displacement  in  the  r9  plane  for 

r<rg,  coupled  to  a rapid  flow  around  the  edge  of  the  cylinder  defined  by 

r~r  . This  motion  is  illustrated  in  Fig.  (VIII  2).  It  is  instructive 
s 

to  compare  Fig.  (VIII  2)  with  Fig.  (VII  1).  Clearly  these  two  types  of 
motion  are  quite  similar;  in  the  former,  a magnetized  fluid  collides  and 
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bounces  off  a stationary  fluid,  while  in  the  latter,  two  fluids  with  equal 
and  opposite  velocity  collide  and  bounce  off  each  other.  Since  the  fluid 
motions  are  so  similar,  we  might  expect  that  this  m = 1 motion  in  a cylindrical 
plasma  is  also  unstable  if  finite  resistivity  (and  therefore  magnetic 
reconnection)  is  allowed.  Indeed,  as  is  apparent  from  Fig.  (VIII  4), 
if  there  is  shear,  the  frozen  in  magnetic  field  lines  have  sharp  turns 
which  provide  strong  restraining  forces.  The  presence  of  resistivity  allows 
the  lines  to  break  and  reconnect,  so  that  they  would  appear  as  in  Fig.  (VIII  6). 
The  field  in  Fig.  (VIII  6)  (analogous  to  that  in  Fig.  (VII  3))  has  fewer 
sharp  turns  and  thus  provides  weaker  restraining  forces.  Thus  for  nonzero 
resistivity,  there  is  more  free  energy  to  drive  the  instability. 

Simpler  diagrams  illustrating  the  same  point  are  shown  in 
Figs.  (VIII  7 and  8).  For  purely  two  dimensional  m = 1 motion  in  a 
cylindrical  plasma  with  = 0,  the  frozen  in  field  lines  in  the  (r0) 
plane  are  shown  at  four  times  in  Fig.  (VIII  7)  (which  is  analogous  to 
Fig.  (VII  2)).  Clearly  sharp  comers  develop  where  the  field  has  to  stretch 
to  follow  the  fluid  motion.  However  if  non-zero  resistivity  is  allowed, 
the  fields  here  will  break  and  reconnect  so  that  the  field  pattern  at 
t = t^  appears  as  in  Fig.  (VIII  8) , (analogous  to  Fig.  VII  (3  )).  In  these 
figures,  the  dots  are  the  magnetic  0 point  singularities  (nulls);  the  two 
dimensional  field  lines  circle  these  points.  For  the  case  of  ideal  MHD, 
no  new  null  point  can  be  produced  and  as  complicated  as  the  field  line 
gets,  it  circles  the  dot  only  once.  However  if  resistivity  is  present, 
the  topological  constraint  is  relaxed  and  new  nulls  (also  called  islands) 
can  be  produced.  In  this  case  an  island  is  produced  near  the  cylinder 
wall  opposite  from  the  direction  of  flow. 
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If  ^ 0,  the  field  surfaces  are  three  dimensional.  Figures 
(VIII  7 and  8)  can  be  regarded  as  projections  of  the  field  lines  back 
into  the  Z = 0 plane  for  ideal  and  resistive  MED.  The  dots  then  correspond 
not  to  nulls,  but  now  to  magnetic  axes.  In  ideal MHD  in  three  dimensions 
the  topological  constraint  preserves  the  number  of  times  one  field  line 
winds  about  another.  However  if  resistivity  is  present,  new  magnetic 
axes  can  be  formed  and  the  field  lines  can  start  to  wind  around  new  axes. 

Now  let  us  qualitatively  discuss  resistive  m = 1 instabilities. 

As  shown  in  Section  (VIII  c) , the  velocity  is  constant  for  r > rg.  Thus, 
for  ideal  MHD,  the  expression  for  has  the  functional  form 


r <r5 

r > Hs  (VIII  34) 


Since  B^  vanishes  at  r = rg,  the  quantity  A from  Chapter  VII  is  infinite 
and  the  theory  of  tearing  modes  described  there  does  not  apply.  The 


growth  rate  must  be  calculated  by  solving  the  coupled  equations  . "3 


c 


C + 


Ifl  K 3r 

Sir  2 r 


(VIII  35) 


subject  to  proper  boundary  conditions. 

We  will  not  follow  this  route,  but  instead  will  derive  a qualitative 
expression  for  the  growth  rate  by  balancing  power  released  in  the  outer 
MHD  region  with  Ohmic  power  dissipated  in  the  inner  region  near  the 
singular  surface.  Since  the  motion  is  two  dimensional  and  incompressible, 
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the  velocity  is  the  curl  of  a vector  potential  Q where 
Q-  -trrVr  tz 


(VIII  36) 


Then,  as  calculated  in  Chapter  VII,  the  power  liberated  per  unit  length  is: 

(VIII  37) 


fF=  \(i  r Vr*)  ■ ^ X f 


A 


The  Z component  of  V x F was  calculated  in  Section  VIII  A and  the  result 
is 

P = -J*  r<Ar  — j c Jr  Vr  L r * *r  r J 


(VIII  38) 


The  quantity  in  the  curly  brackets  vanishes  in  the  outer  MHD  region, 

exactly  as  in  Chapter  VII.  However  near  the  singularity,  the  bracket  is 

no  longer  zero,  since  ideal  MHD  with  y * 0 is  no  longer  an  accurate 

description  of  the  plasma.  In  a narrow  region  of  width  xc  the  slope  of 

d^Br 

goes  from  zero  to  some  constant  value.  Thus  2 is  the  dominant 
term  in  the  bracket.  The  functional  ibrm  of  B^(r)  is  shown  in  Fig.  (VIII  9) 


as  the  solid  curve,  while  the  ideal  MHD  solution  is  shown  as  the  dotted 

d^B 

curve.  Clearly,  B^  and  have  the  same  sign,  so  that  Pp  is  positive 


and  energy  is  available  to  drive  instability.  If  the  value  of  B^  at 
r » rQ  - x,  is  denoted  B^,  Eq.  (VIII  38)  yields  the  approximate  expression^ 


?F  ~ 


» r5J  | Br< 


(VIII  39) 


The  Ohmic  power  per  unit  length  dissipated  in  the  singular  layer  is  given 
roughly  by 


(VIII  40) 
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where  J z is  the  current  flowing  in  the  singular  layer. 


However 
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(VIII  41) 


where  the  first  relation  in  Eq.  (VIII  41)  comes  from  Maxwell's  current 
equation  and  the  second  from  V*B  * 0.  Then  combining  Eq.  (VIII  39  and  40) 
we  find 


Virxt 

The  only  problem  now  is  to  determine  x in  terms  of  y and  n.  The 


(VIII  42) 


procedure  here  is  exactly  the  same  as  in  Chapter  VII.  Namely,  substitute 
d2B 

for  the  — jpf-  term  in  Eq.  (VIII  35a)  from  Eq.  (VIII  35b).  Then,  exactly 


as  in  Eq.  (VII  23), 

A ~ 


p *rc' 

rT 


* 


(VIII  43) 


Thus,  one  can  derive  in  a fairly  simple  way  the  basic  results  of  m = 1 
tearing  modes. 


Ill 


VIII  F - Double  Tearing  Modes 


We  conclude  this  chapter  with  a study  of  double  tearing  modes; 

that  is  tearing  modes  which  are  excited  when  there  are  two  zeros  of  F 

near  each  other.  For  tokamak  ordering,  this  generally  means  a current 

density  which  is  not  a monatonically  decreasing  function  of  radius,  but 

which  peaks  at  some  radius,  in  other  words  a skin  effect.  This  can  be 

seen  from  Fig.  10  a and  b) . The  solid  curves  are  normal  current  density 

and  poloidal  field  as  a function  of  radius.  The  dotted  curves  show 

the  effect  of  a narrow  skin  layer  on  toth  J and  B„.  Since  F = B„  - kBz, 

9 >r 

F clearly  has  basically  the  same  behavior  near  rc  as  the  dotted  curve 
in  Fig.  (VIII  10b). 


To  investigate  the  stability  of  such  a configuration,  first  look 
at  the  regions  of  the  plasma  where  ideal  MHD  is  valid.  Using  the  definition 
of  F and  Maxwells  current  equation,  Eq.  (VIII  28)  can  be  reduced  to 

_F  _ _Br__^r3—  * 


r1  Jr  Jr 


*•  A A 


7rf  - / „ 

F 


(VIII  46) 


This  equation  is  to  be  solved  subject  to  appropriate  boundary  conditions 
at  r * 0 and  r * a. 

If  r and  r are  the  positions  where  F vanish,  and  if 
cl  c2 

r - r <<  r , then  r 'v-  r , the  second  term  on  the  right  hand  side  of 
Cl  C2  c c* 

Eq.  (VIII  46)  is  much  smaller  than  the  first,  and  can  be  neglected.  Making 
this  approximation,  one  can  easily  construct  a solution  to  Eq.  (VIII  46) 
which  satisfies  boundary  conditions  at  r * 0 and  r = a and  is  everywhere 


continuous.  The  solution  is 


r<  rC3 

rc  a < r < rc, 
rc,<  r 


(VIII  47) 
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This  is  analogous  to  the  solution  for  for  m = 1 tearing  modes. 

If  p , n and  F"  have  the  same  value  at  r = r and  r , an  analysis  like 

o ci  C2 

that  in  the  previous  section  gives  Eq.  (VIII  44  and  48)  for  growth  rate  and 

size  of  singular  region.  (Of  course  x^  <<  rcl-rc2  is  also  assumed).  Notice 

that  the  growth  rate  for  a skin  current  goes  as  n1/3,  whereas  for  a normal 

3/5 

tearing  mode  it  goes  as  n . Since  n is  a small  quantity,  a double  tearing 

mode  grows  much  faster  than  a conventional  tearing  mode. 

Since  the  ideal  MHD  relation  between  V and  B is  yB  = iFV  , the 

r r r r 

solution  for  B^  in  Eq.  (VIII  97)  gives 


r<r. 


Vr-  < co*s rci<  r<  rc, 

Jo  rC(  < r 


(VIII  US) 


The  flow  pattern  is  then  that  which  is  characteristic  of  tearing  modes, 
namely  a fluid  collides  with  and  recoils  off  a stationary  fluid.  Such 
a flow  pattern  is  illustrated  in  Fig.  (VIII  11 ^ for  the  case  of  m * 4. 


113 


Fig.  VIII  2 — The  flow  pattern  at  four  axial  positions.  The  dots  labeled  1 , 2,  and  3 are 
intersections  of  a field  line,  having  the  same  pitch  as  the  flow,  with  these  four  planes  at 
three  different  times.  The  points  labeled  A are  intersections  of  a field  line  having  q < 1 
with  these  four  planes.  At  Z = 0,  point  A is  closer  along  a stream  line  to  the  separation 
region.  Therefore  it  will  enter  this  region  first  and  will  race  around  the  cylinder  before 
it  does  so  at  other  axial  positions.  This  will  lead  to  considerable  bending  and  stretching 
of  the  field  line. 
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Fig.  VIII  3 — A three  dimensional  view  of  the  motion  of  two  neighboring  field  lines 
having  q = 1 (same  pitch  as  the  flow).  The  field  lines  are  red  and  the  flow  pattern  is 
blue.  A solid  field  line  means  it  is  near  the  front  of  the  cylinder,  dotted  means  near 
the  back.  At  t = 0 (a)  the  field  lines  in  front  of  the  cylinder  are  converted  towards 
the  separation  region.  At  tj,  (b)  they  both  enter  the  separation  region  and  whip  around 
the  cylinder  in  opposite  directions.  At  t2  (c)  they  are  continuing  around  the  edge. 
Finally  at  t3  (d)  they  have  both  re-entered  the  main  rigid  body  flow,  and  now  are  in 
back  of  the  cylinder. 
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(a)  (I))  (c)  (cl) 


Fin.  VII M A throe  dimensional  plot  of  a field  line  having  higher  pitch  (q  < 1 ) than 
the  flow  pattern.  At  t = 0,  the  field  line  wraps  around  the  inner  cylinder  shown  in  (a). 
Alt  - t] , (h)  this  cylinder  convents  like  a rigid  body  (in  the  rO  plane)  with  the  main 
flow  and  the  field  line  continues  to  wrap  around  it.  At  t = t2  (c),  the  field  line  enters 
the  separation  region  at  one  point  and  begins  to  whip  around,  while  most  of  the  field 
line  just  continues  to  move  as  a rigid  body.  At  t = 13,  (d)  more  of  the  field  line  has 
entered  the  separation  region  and  the  first  part  of  the  field  line  to  enter  the  separation 
region  has  re-emerged  out  the  other  side  of  the  cylinder.  Clearly  the  field  line  is  con- 
torted much  more  than  in  Fig.  VIII  3. 


Fig.  VIII  7 — A two  dimensional  field  line  in  the  r,d  plane  which  is  frozen  into  the 
flow  given  by  Eq.  VIII  25.  The  field  line  is  shown  at  four  different  times.  Notice 
that  it  is  rapidly  converted  around  the  separation  layer  and  re-emerges  into  the  main 
flow  on  the  other  side  of  the  circle.  At  t = t3  fields  of  opposite  sign  are  forced  next 
to  each  other  in  the  separation  layer  and  along  a horizontal  diameter  on  the  left  side 
of  the  circle. 


Fig.  VIII  9 — A plot  of  Br  as  a function  of  r for  an  m 
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IX.  Instabilities  in  a Toroidal  Plasma 

In  this  chapter,  we  just  touch  upon  the  question  of  MHD  instabilities 
in  toroidal  geometry.  The  actual  calculations  of  relevant  instabilities  has 
many  more  mathematical  complications  than  the  equivalent  calculations  in 
slab  or  cylindrical  geometry.  The  reason  of  course  is  that  toroidal  plasmas 
are  inherently  two  dimensional.  For  instance  if  the  toroidal  co-ordinates 
are  as  shown  in  Fig.  IX.  1,  the  equilibrium  is  symmetric  in  toroidal  angle  l 

S,  but  depends  on  the  two  variables  r and  9.  Thus  calculating  the  stability 
involves  a two  dimensional  eigenvalue  equation  in  r and  9,  rather  than  one 
dimensional  calculations  (in  r for  cylindrical  geometry,  or  in  x for  slab 
geometry)  as  done  in  the  preceeding  chapters. 

We  will  not  attempt  such  a two  dimensional  calculation  here;  rather 
will  examine  first  how  modes  in  cylindrical  plasmas  are  affected  by  toroidicity, 
and  second  will  show  that  there  are  instabilities  in  a toroidal  plasma  which 
have  no  analog  in  cylindrical  geometry.  In  tokamak  ordering,  the  results 
are  that  free  surface  modes  and  tearing  modes  exist  as  in  cylindrical 
geometry.  The  pressure  driven  modes,  described  in  Chapter  VIII  b,  are 
stabilized  of  q > 1.  Finally  a new  mode,  usually  called  a balooning  mode, 
can  exist  in  the  regions  of  locally  bad  curvature.  These  balooning  modes 
give  rise  to  a pressure  limit  for  tokamak  plasmas. 


We  now  examine  these  points  more  thoroughly.  In  the  tokamak  ordering, 

a<<R,  B <<B  , and  all  toroidal  effects  will  be  small  by  a factor  of  at  least 

o z J 

a/  . However  in  deriving  the  equations  for  modes  with  two  dimensional 
structure  in  cylindrical  geometry  (Eqs.  (IV  35”  ) or  (VIII  'J  )),  all 
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terms  of  order  a,  were  neglected.  Therefore,  any  instabilities  which 
R 

exist  in  this  approximation  for  a cylindrical  plasma  should  also  exist  in 
toroidal  geometry  if  a<<R.  The  relevant  instabilities  are  then  the  free 
surface  modes  described  in  Chapter  IV  and  the  tearing  modes  described  in 
Chapter  (VIII.  d-f).  However,  as  pointed  out  in  Chapter  (VIII.  a),  a 
plasma  with  no  free  surface  and  with  zero  resistivity  is  stable  to  modes 

with  two  dimensional  structure.  Now  let  ib  examine  the  pressure  driven  t 

modes  discussed  in  Chapter  (VIII.  b) . These  modes  were  driven  by  a pressure 

gradient,  and  as  is  apparent  from  Eq.  (VIII. /«?c) , the  driving  term  is 

multiplied  by  k2r2  ^(^)2.  However  toroidal  effects  are  also  of  this 

order,  so  that  these  pressure  driven  modes  can  be  strongly  affected. 

Calculation  of  the  exact  stability  condition  is  mathematically  involved  and 
was  first  given  by  Mercier  and  also  by  Shafranoy^ 

We  will  give  here  a plausible  argument  tfiich  gives  the  correct  result,  namely 

rB 

zo 

the  pressure  driven  modes  in  a tokamak  are  stabilized  if  q = — — >1. 

KB. 

9o 

Recall  that  Chapter  V speculated  that  the  effective  gravity  arose  from 
forces  acting  on  the  fluid  which  originate  in  the  geometrical  complexities 
of  Bq  which  had  no  analog  in  slab  geometry.  This  was  confirmed  in 
Chapter  VIII.  b where  the  gravity  was  related  to  the  curvature  of  field 
lines  in  cylindrical  geometry.  We  will  use  this  result  to  calculate  the 
appropriate  gravity  for  a toroidal  plasma.  Here,  there  are  three  geometrical 
effects  which  give  rise  to  an  average  force  having  no  analog  in  slab  geometry. 

The  first  and  second  are  forces  due  to  the  curvature  of  both  the  toroidal 
and  poloidal  field,  the  third  is  the  outward  force  due  to  the  fall  off  in 
toroidal  field  strength  with  major  radius.  This  latter  is  clearly  a 
geometrical  effect  since  it  is  not  balanced  directly  by  a change  in  plasma 
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pressure.  The  average  force  on  a particle  arising  from  the  curvature  of 
the  field  lines  is 

Fc= 


(8.- 1 ) ?. 
8.‘ 


(IX.  1) 


while  the  force  on  a particle  resulting  from  the  gradient  of  the  strength 
of  the  toroidal  field  is 

r s M V I Br.  I 

* * 'TbTT  (ix*  2) 

Adding  up  the  forces  on  all  particles,  and  using  the  fact  that  ("^V^2)-  =■ 

< V2  > = V2  for  an  isotropic  plasma,  we  find  that  the  effective  gravitational 

it  i 

force  g is  given  by 

- _ v,1  ((*■*) & + vl&J 

^ ‘ { mi1  1 3J 

Let  us  first  calculate  the  (B  *V)B  term  above.  To  do  so,  let 

-o  - -o 

introduce  cylindrical  co-ordinates,  shown  in  Fig.  (IX.  2).  Here  f is  the 
same  angle  as  in  Fig.  IX.  1.  The  components  of  the  field  in  cylindrical 
co-ordinates  are 

1? 
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(IX.  3) 


R - B — 

' D'°  R 


(R-R.) 

Be“  />-/?.)'  + z1] 


V 


= 3. 


B°  [(H-  1T.)V  *‘J 


00 


(O 


(IX.  4) 
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Then  after  some  algebra,  one  can  make  use  of  the  expression  for  (Bq*V)B 
in  cylindrical  co-ordinates  to  give  the  result 


o 


- Be.  2 

(R-  Z1 


(IX.  5) 


where  we  have  included  only  the  components  in  the  Rz  plane.  In  Eq.  (IX.  5), 

the  first  term  multiplying  i and  the  i part  clearly  arise  from  the  curvature 

“K  “2 

of  the  poloidal  field;  the  second  term  multiplying  i arises  from  the 

— K 

curvature  of  the  toroidal  field. 

Now  pressure  driven  modes  can  occur  when  the  density  (or  pressure) 
gradient  is  anti  parallel  to  g.  However  the  density  gradient  is  in  the 
negative  r direction  (here  r is  given  in  Fig.  IX  1 ).  Using  the  fact 


that 
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(IX.  7) 
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Similarly 
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Thus  our  final  expression  for  g is 
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(IX.  9) 


Assuming  that  B.  <<B  , we  find 

9o  20 


Coo-  8 


K + rcosl 


(IX.  10) 


The  first  term  in  the  bracket  of  Eq.  (IX.  10)  is  simply  from  the  curvature 
of  the  field  line  in  cylindrical  geometry.  The  second  is  from  the  curvature 
and  gradient  of  the  toroidal  field.  It  is  a larger  term  than  the  first,  but 
it  nearly  averages  to  zero  over  poloidal  angle  9.  In  order  to  see  the  effect 
of  this  term,  expand  the  denominator  in  powers  of  rc°s^  and  average  over  9. 
The  result  is 


< 3-i->  = 


(IX.  11) 


dn  dP 

Thus  for  q>l,  the  averaged  g changes  sign  and  is  negative.  Since  — — and  — - 


dr 


dr 


are  also  negative  the  effect  of  the  gravity  is  no  longer  destabilizing. 
Therefore,  tokamaks,  which  are  characterized  by  q>l,  are  stable  to  pressure 
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driven  modes.  Reversed  field  pinches  however,  which  have  q<l,  are 
potentially  unstable  to  pressure  driven  modes  in  ideal  tlHD  if  the  shear 
is  not  great  enough,  and  are  also  potentially  unstable  to  resistive  g 
modes  no  matter  what  the  shear  is. 

We  close  this  chapter  by  qualitatively  examining  the  problem  of 
balooning  modes  in  tokamaks.  Notice  that  the  gravity,  g'i^  in  Eq.  (IX.  10) 

is  the  sum  of  two  terms,  one  of  which  is  constant  along  the  field  line  and  ( 

one  which  varies  along  a field  line.  The  latter  nearly  averages  to  zero 
along  a field  line  because  of  the  cos  9 in  the  numerator.  However,  in  tokamak 
ordering  the  maximum  value  of  this  latter  term  is  much  larger  than  the  value 

2Rq2 

of  the  former  (*~  ^ ■ larger,  to  be  precise).  This  naturally  brings  up  the 
following  question:  Can  a mode  be  localized  along  a field  line  so  that  it 
only  sees  the  positive,  destabilizing,  values  of  g*i  (around  9~0)  and  not 
the  negative,  stabilizing,  values?  We  will  see  that  such  instabilities  can 
indeed  be  generated  if  the  plasma  pressure  is  too  high.  These  are  usually 
called  balooning  modes,  probably  in  comparison  to  the  behavior  of  a baloon 
of  varying  thickness.  Imagine  a baloon  whose  thickness  varies  around  its 
surface,  also  imagine  a second  baloon  of  uniform  thickness  equal  to  the 
average  thickness  of  the  first.  Say  that  someone  is  too  weak  to  pop  the 
second  baloon;  he  might  still  be  able  to  How  up  and  pop  the  first,  because 
it  blows  out  much  more  easily  where  the  surface  is  thin. 

This  is  then  analogous  to  an  MHD  instability  in  a tokamak  picking  out  - 

the  most  unstable  (outer)  region  and  localizing  itself  there.  Of  course, 
it  is  not  so  simple  for  a mode  to  localize  itself  along  a field  line.  Any 
disturbance  which  is  localized  along  a field  line  is  also  setting  up  shear 
Alfven  waves,  which  tend  to  give  a positive  frequency  vj^L2  where  L is  the 
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characteristic  length  of  the  disturbance  along  the  field  line  and  V,  is  the 

A 

Alfven  speed.  The  distance  L is  given  roughly  by  (2ir ) — 1 times  the  distance 
along  the  field  line  between  0 = - /£  and  9 3 or  L a Rq^-  Ttie  maximum 
growth  rate  is  given  roughly  by  the  maximum  value  of  g divided  by  the  gradient 
scale  length,  taken  here  as  / The  plasma  can  only  be  stable  if  the 
outward  motion  forced  by  the  gravity  does  not  overcome  the  resistance  to 
field  line  bending,  or  if 


This  gives  rise  to  an  approximate  limit  on  the  total  8 of  a tokamak  of 
V2 

3 = II  < 2l- 
® V2  - Rq2  ’ 

Now  let  us  show  how  this  basic  result  can  be  obtained  more  rigorously. 
If  we  assume  slab  geometry  with  all  perturbed  quantities  varying  as 
(z)exp  i ky,  Bq  uniform  and  in  the  z direction,  the  density  gradient 
in  the  x direction,  and  g varying  with  z,  an  analysis  similar  to  that 
which  led  to  Eq.  (V.  n ) gives  the  result 


(-f 


fW  2 z‘ 


-4 


= 0 

x z * J 


313  V2  cos  0 

Specializing  to  a tokamak,  where  — = — — and  g(z)  = 2 

d Z Kq  o\j  K 

that  Eq.  (IX.  12)  reduces  to 

l 


(IX.  12) 


we  f ind 


~ 7 t * j 0 

*/  M 


K U 


V ,ie)  = o 


(IX.  13) 


where  Ln  is  the  gradient  scale  length. 
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Notice  that  the  last  term  in  the  parentheses  is  periodic  in  9 with 


period  2tt.  However  the  solution  for  V (9)  must  also  be  periodic  with  the 

same  periodicity,  since  Vx  must  be  a single  valued  function  of  0.  Thus  y 

is  determined  by  the  condition  that  solutions  to  Eq.  (IX.  13)  exist  which 

are  periodic  in  9 with  period  2tt.  An  approximate  criterion  for  unstable 

roots  can  be  derived  as  follows.  If  y = 0,  the  solutions  to  Eq.  IX  are 

oscilliatory  between  - ^ < 9 < and  are  exponentially  growing  and 

damping  for  <|9|<ff*  Imagine  a solution  for  V localized  between  - ^ 

and  and  which  damps  to  a very  small  value  between  say  ^<d<-n.  Then 

the  solution  for  between  *•  and  does  not  affect  the  solution  for 

3tt^  5 it, 

Vx  betweem  say  ^ and  The  condition  for  a solution  to  Eq.  (IX.  4) 

with  zero  y is  simply  the  condition  for  localized  roots  between  the  turning 


points,  or  according  to  WKB  theory, 


I 


JK 

Ln 


'/a. 


Cov* 


= 1? 
/a. 


(IX.  13) 


The  eigenfunction  is  as  shown  in  Fig.  (IX.  3).  Notice  that  between 

< 9 < it,  gets  so  small  that  it  does  not  affect  the  value  of  Vx  in 


JJ  ' 5lL 


< 9 < 3^  » where  Vx  gets  large  again.  Using  the  fact  that  J*  ^cos  d6  - 


we  find  that  the  plasma  is  unstable  if 

v; 


(IX.  16) 


This  condition  then  puts  a limit  on  the  maximum  pressure  which  can  be 
stably  confined  in  a tokamak.  If  the  plasma  is  unstable,  growth  rates 
can  be  calculated  in  an  analogous  way  by  setting 
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The  problem  of  balooning  modes  in  sheared  fields,  and  with  nonzero 
resistivity  is  much  more  complicated  and  is  still  under  active  investigation. 
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Fig.  IX  2 — The  toroidal  co-ordinate  R,  z and  f. 


Fig.  IX  3 — The  eigenfunction  (dotted  curve)  and  local  wave  number  squared 
(solid  curve)  for  ballooning  modes  in  a tokamak. 
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